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I 

INTRODUCTION 

The feature of mathematics which we may call scientific approxi- 
mation did not receive the recognition among the ancient Greek 
geometers that has been accorded it by mathematicians of modern 
times. Since their work in pure mathematics was based on geometry, 
the Greek scholars evaluated irrational roots as readily as rational 
roots by the construction of lines. But when it came to applied 
mathematics, nvunerical approximation was often not only desirable 
but absolutely necessary; and so we have some remarkably close 
approximations by Archimedes (c. 225 B. c.) in his work in mechan- 
ics, by Heron of Alexandria (c. 200 A. D.) in his work as surveyor, 
and by Theon of Alexandria (c. 375 a. d.) in his astronomical com- 
putations. Archimedes even used numerical approximations in the 
domain of pure geometry in his Measurement of the Circle. On the 
whole, however, the Greek mind did not take to this form of mathe- 
matical activity. 

The Hindus and the Arabs used approximation methods more 
extensively than the Greeks, and notable work was done by Brah- 
magupta (c. 510 A. D.), Bhaskara (c. 1150 A. d.), and an unknown 
Arab computer mentioned by Chelebt. However, it was during the 
Renaissance, when the newly discovered general formula for solv- 
ing the cubic was found to be inoperative for the "irreducible case," 
that approximation methods became a vital problem, and for two 
centuries the keenest mathematical minds worked to find smooth, 
effective methods for approximating cube and higher roots. 

The purpose of this research is to trace the history of the different 
methods of approximating roots of numerical higher equations that 
were used up to 1819, the date of the publication of Horner's method; 
we purpose to trace their early beginnings in finding the roots of 
numbers and in solving incomplete equations, and to watch their 
growth into systematically developed general methods for solving 
complete equations of any degree. Pure trial and error methods, 
like some of the methods of the Egyptians and, in modern times, 
that of Junge, will not be taken up, nor shall we discuss such 
ephemeral methods as those of Fontaine, Collins, and some of 



2 Approximating Roots of Numerical Higher Equations 

DeLagny's. We shall, moreover, limit the investigation to algebraic 
methods in algebraic equations of one unknown, and only inci- 
dentally shall we touch upon transcendental equations and geo- 
metric and trigonometic methods. 

In searching for literature on this subject the author found only 
two articles, outside of the encyclopedias, in the nature of surveys: 
one by Augustus De Morgan in the Companion to the British Almanac 
(1839), entitled "Notices of the Progress of the Problem of Evolu- 
tion" and dealing with the extension of the Hindu method of evo- 
lution to the solution of numerical equations by Vieta and Horner; 
and a more extended discussion by F. Cajori in an article entitled 
"A History of the Arithmetical Methods of Approximation to the 
Roots of Numerical Equations of one Unknown Quantity," in the 
Colorado College Publication (Colorado Springs, Colorado) for 1910. 
In the present survey we purpose to study the origin and growth of 
the various methods by bringing out the illustrations and the details 
of explanation as far as possible from the original sources. 

The writer has had unusual opportunities for studying such source 
materials in four notable libraries, — the Columbia University 
Library, the Teachers College Library, the New York City Public 
Library, and the private library of Professor David Eugene Smith, 
with its unique collection of editions of the Renaissance and early 
post-Renaissance periods. By this means he has had access to the 
early editions and standard translations of Euclid, Archimedes, 
Heron, Theon, Brahmagupta, Mahavira, Bhaskara, Omar Khay- 
yam, as well as the latest brochures and translations of Chinese and 
Japanese writers on mathematics. He has also been able to consult 
such works as those of Leonardo of Pisa, Chuquet, Pacioli, Bom- 
belli. Cardan, Stevin, Pitiscus, Vieta, Oughtred, Hume, Herigone, 
Wallis, Newton, Raphson, Colson, Rolle, DeLagny, Halley, Taylor, 
Waring, Euler, Lagrange, and many others who wrote upon the 
subject previous to Horner. 



II 

PRE-ALGORITHMIC METHODS 

I . EGYPTIAN AND BABYLONIAN WORK IN FINDING ROOTS 

We formerly thought that the root concept originated among 
the Greeks, but recent discoveries show us that the finding of roots 
of numbers and the solution of quadratic equations, and even the 
approximation methods employed in this work, go as far back as 
the early Egyptian and Babylonian civilizations. 

On the Senkereh tablets (c. 2000 b. c.) are tables of squares and 
cubes, which shows that the Babylonians had at least an indirect 
notion of square and cube roots. ^ Three quadratic equations are 
known to have been studied in this period. The first one was made 
known when Griffith (1897) published the mathematical papyrus 
found by Petrie in Kahun.* It deals with areas and requires the 
solution of the equations xy = 12 and x :y = i i^, stated, of 
course, in rhetorical form. In 1900 Schach discovered in a Berlin 
papyrus a second problem in quadratics, requiring the solution of 
the equations x^ + y^ = 100 and x -.y = i : ^. The third equation 
was found in the Kahun papyrus by Schach in 1903. It requires the 
solution of the equations x^ + y^ = 400 and x :y = 2 : i}4. The 
ancient mathematician solved it by letting x = 2, y = i}4; this 
gives x^ + / = 6><; since V6K = 2>^ and 2>^ = Ys oi 20, he 
found that a: = 2-8 = 16 and y = i>^-8 = 12. 

The same method, commonly called the method of False Position, 
is also used in solving the other two equations. 

No irrational roots occur in these equations. But there are indi- 
cations ' that the Egyptians had a definite way of approaching the 
square root of non-square nimibers. However, they seem to have 
been unaware of their irrational quality. 

2. EARLY CONCEPTS OF HIGHER ROOTS AND IRRATIONALS 
AMONG THE GREEKS 

That some roots are irrational was first recognized by the Greeks. 
From measuring the areas of commensurable squares and rectangles 

1 M. Cantor, Vorlesungen iiber Geschichte der Mathematik, Vol. I (hereafter desig- 
nated Cantor, I), Leipzig, 1907. P- 28. 

' M. Simon, Geschichte der Mathematik im Altertum, Berlin, 1909, pp. 41-42. 
= Simon, Geschichte, pp. 43-53- 
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in terms of unit squares they proceeded to find the areas of non- 
rectangular and irregular figures, such as the triangle and the circle. 
This they did by finding an equivalent square (Quadrature). Since 
they neither possessed algebraic symbols nor used the equational 
form of statement, they found the unknown side by geometrical 
construction. Some of the sides were discovered to be incommensu- 
rable quantities. As early as 525 B. c. Pythagoras recognized the 
diagonal length of the unit square to be an irrational number. 
Theodorus of Cyrene (c. 410 b. c.) and Plato (c. 380 b. c.) discussed 
the irrational nature of the expressions Vs, Vs, y/j, ■ ■ ■ "s/17. 
In Book X of his Elements Euclid (c. 300 B. c.) gave a comprehensive 
exposition of the theory of irrationals. Up to the time of Archi- 
medes square roots were almost universally evaluated by geometri- 
cal constructions; by him they were often evaluated numerically 
by a process of approximation. 

The Greeks seem also to have been the first people to conceive 
directly the notion of cube root.'' This also arose from geometrical 
problems, and, like the square root, was approached geometrically. 
According to the familiar story the god's demand of an altar twice 
as large as the old one necessitated the solution of the equation 
X' = 2a'. This was probably solved by Menaechmus (c. 365 B. c.) 
by the aid of conic sections. Archimedes' famous problem of the 
section of the sphere led to an equation of the type x^ — ax^ -|- J = o, 
which was also handled geometrically. Heron of Alexandria (c. 
200 A. D.) was the first one to make a practise of evaluating cube 
roots numerically; this he did in the course of his practical work by 
a process of approximation. 

3. THE APPROXIMATIONS OF ARCHIMEDES AND HERON OF ALEXANDRIA 

In his Measurement of the Circle Archimedes (287-212 B. c.) 
gives among others the following approximations:*^ 
1351 -^ /T ^ 265 



^„ > \/3 > ; V349450 > 591-; V5472132T5 > 2339-- 

780 153 8 4 

Nowhere has he disclosed his process, and many conjectures and 
many ingenious restoration formulas have been advanced as ex- 

' Tropfke, Geschichte der EUmenlar-Mathematik I (i), Leipzig, 1902, pp. 158-61; 
207-10; 223-24. 

' Cantor, I {1907), p. 316. 
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planations.* We have one bit of information. In discussing Heron 
of Alexandria (c. 200 a. d.), who Hved four centuries later, Eutocius 
tells us that he found square and cube roots by the same method as 
Archimedes. That remark, together with testing for certain recon- 
struction formulas, seems to indicate that at times he used methods 
of Double False Position,' as in evaluating the expression V2^ and 

at other times he used the formula VN = - ( a + — ),asin find- 

2 \ a / 

ing the square root of 3. 

Now, we have Heron's processes described in his own words. A 
tenth century manuscript of Heron's Metrica was recently discov- 
ered by Dr. R. Schone in the Serailo Library in Constantinople. In 
describing how to find the square root of 720, the area of a certain 
triangle. Heron says: 

"Now, since 720 has no rational square root, we shall find the 
root differing from it by a very small error. Since the next larger 
square to 720 is 729 and its side is 27, we divide 720 by 27. The result 
is 26%. To this add 27 ; that gives 53^^. One-half of this is 26 }4 H- 
Therefore the next root of 720 is 26 }4 K- Multiplying 26 }4 }i 
by itself gives 720V86, so that the error amounts to only Vm of unity. 
Should we wish, however, to obtain an error smaller than Vse, we 
may put in place of 729, the value we have just found, 720 Vse, and 
if we do this, we shall find an error much smaller than Vse." 

For want of a name we shall call this the "method of averages." 
If a is a rational root of N, then N-i-a = a. But if it is an approxi- 
mation, and a little too large, as in Heron's example, then 
N -i- a = fli, a little smaller than the exact root. Heron takes as 

the final approximation -(a + ai) = -la + — )• 

2 2\ a/ 

For approximating the cube root * Heron used the method of 

Double False Position. He explains this process as follows: "We 

shall now explain how to find the cube root of 100 units. Take the 

two cubes nearest to 100, one larger, the other smaller. Note how 

'Cantor, I (1907), pp. 315-17; consult S. Giinther, "Quadratische Irrationaliteten 
der Alten," in Abhandl. zur Geschichte der Math., Vol. IV. 

' G. Wertheim, "Heron's Ausziehung der irratlonalen Kubikwurzeln," in Zeitschrift 
fiir Math. u. Phys., XLIX (1899), Hist. Lit. Abt., pp. 1-3. 

' Curtze, in Zeitschrift fiir Math. u. Phys., XLII (1897), p. 117-119; Wertheim, in 
XLIV (1899), pp. 1-3; Bibliolheca Malhemalica (hereafter called Bibl. Math.), third 
series. VIII (1907-1908), p. 412; Cantor, I (1907), p. 374. 
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much larger the first one is, namely 25, and how much smaller the 
second one is, namely 36. Then multiply 36 by 5; the result is 180. 
To this add 100, which gives 280 [and divide 180 by 280] ; the result 
is %4. Add this to the root of the smaller cube, namely 4; the result 
is 4 %4. Which is the value of the cube root of 100 units as nearly 
exact as possible." 

The ambiguity of the nimibers in this example (5 maybe -^^125 
and also ^125 - 100, etc.) excludes a precise knowledge of Heron's 
process; but it seems to differ from those used by the medieval 
writers. Curtze, Wertheim, Cantor, and Enestrom have all worked 
out different reconstruction formulas for Heron's process for 
approximating cube roots. 

4. theon's method of exhaustion 
The first instance in extant writing of a method of finding square 
root similar to our present method is found in the writings of Theon 
of Alexandria (c. 365 a. d.). He gives a very full description of 
how to evaluate irrational roots in sexagesimal fractions.' It is a 
geometrical process referred to a figure. In it he applies a process 
of exhaustion to Euclid's geometrical counterpart of the expression 
(a + xY = a^ + 2ax + x^. This is followed by a generalizing rule. 
Expressed in modern analytical form, if a^ -\- b = {a -\- xY = 

a^ + 2ax + x^ and a is a first approximation, then x = — approxi- 

2a 



mately, and Va^ -\- b = a -\ , approximately. In his illustrative 

2a 

example he finds the square root of 4500° to be 67° 4' 45". 

Diophantus did brilliant work in approximation, but this was 

chiefly in connection with inequalities and indeterminate equations, 

both of which are outside of our topic. 

5. SUMMARY 

1. We find the first concept of roots in Egypt or in Babylon; in 
Egypt we also find the solution of quadratic equations by single 
false position. 

2. The concepts of irrationality and of higher roots are first found 
among the Greeks. Theirs are also the earliest records of a scien- 
tific notion of limit and of scientific approximation. 

' Gow, History of Greek Mathernatics, Cambridge, 1884, p. 55. 
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3. We find among them the first recorded use of three approxi- 
mation methods: namely, the methods of averages and of double 
false position, by Archimedes and Heron; and the method of ex- 
haustion, by Theon. Theon's method differs from Heron's methods 
in giving sharper relief to the limit idea ; the method of averages and 
the method of double false position really have two variable limits. 
Though differing in principle and operation, Heron's method of 
averages and Theon's method of exhaustion give the same results 
at each step. For, if N = a' + b, then will y/N = Va^ + b = 

2\ a f 2\ a ' 2a 

4. The existence of a root, and never more than one root, was 
taken for granted, among both Egyptians and Greeks. 



Ill 

ALGORITHMIC METHODS FOR APPROXIMATING 
THE ROOTS OF PURE POWERS 

I . THE HINDU METHOD OF EXHAUSTION 

The superior notation and convenient symbolism of the Hindus 
enabled them to use processes of inversion to a degree impossible 
to Theon. In fact that became the favorite line of attack in much 
of their mathematics. Square and cube roots were found by 
Aryabhafa ' {b. 476 A. d.) from the inverse of the formulas 
a" + 2ab + b^ and a^ + sa'b + sab^ + b\ 

When the Hindu decimal notation with its place values came 
into use, the work of finding roots was so arranged that an applica- 
tion of the inversion process would evolve the digits of the root in 
order (whence "evolution" in texts on arithmetic and algebra). To 
effect this the number was divided into periods and the work 
arranged in columns and lines. This "Hindu method," as it was 
often called in the Middle Ages, is described for us by Brah- 
magupta (&. 598 A. d.).^ Sridhara,' in his Triiatikd (1040 a. d.) 
gives the rules in full, as also does Bhaskara (&. 11 14 A. d.).* The 
following is a free translation from Bhaskara, by Taylor: 

"The first place on the right is called ghana or cube ; the two next 
places aghana or not-cube. — Subtract the cube contained in the 
final period from the said period; put down the root of the cube 
in a separate line, and after multiplying its square by three, divide 
the antecedent figure by the result, and write down the quotient 
in the separate line: Then multiply the square of the quotient by 
the preceding number in that line and by three, and after sub- 
tracting the product from the next antecedent figure cube the said 
quotient, and subtract the result from the next antecedent figure. 
Thus repeat the process through all the figures. The separate line 
contains the Cube Root." 

1 Cantor, I (1907), p. 625. 

2 H. T. Colebrooke's translations of Brahmagupta and Bhaskara, London, 1817, 
pp. 279-81. Hereafter referred to as Colebrooke. 

' Kaye, on Trislika, in Bibl. Math. (3) XIII (1912-13), p. 209. 
' Taylor's translation of Lilawati, Bombay, i8i6, p. 20. 
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The Hindu method was adopted by the Arabs, and Omar Khay- 
yam (c. 1 100 A. D.) wrote an exposition of it in a work that has been 
lost. From the Arabs it came to Christian Europe * through the 
writings of Leonardo of Pisa (1202); and the arrangement of the 
work as given in Sacrobosco's Algorismus (c. 1240) remained in use 
until it was supplanted by that of Peurbach in the fifteenth century. 

Our modern arrangement of the work in using the Hindu method 
begins with Peurbach (1423-1461), whose prestige as an arithmeti- 
cian and an astronomer caused his schema to be adopted throughout 
Europe." Slight modifications were added by Chuquet (1484), 
LaRoche (1520), and Gemma Frisius (1540). In the Oriental coun- 
tries an abacus-like arrangement obtained. 



2. FORMULAS FOR IRRATIONAL ROOTS IN THE MIDDLE AGES 

The Hindu view of irrationals differed somewhat from that of the 
Greeks. To the latter, irrational numbers represented incom- 
mensurable lines, and hence they were slow to evaluate them 
arithmetically. To the Hindu an irrational number was "one, the 
root of which is required, but cannot be found without residue."' 
As far back as before 200 A. d., in the Sulvasutras,' we find re- 
markably close approximations for y/2, Vs, and other surds, by 
the altar computers. Thus Baudhayana states that \/2 = 

I -\- - -\- — — . Their method is not known. It does not 

3 3-4 3-4-34 
seem to have been a method of inversion. But later writers ^ used 
the principle of inversion for irrationals almost entirely. 

Among the Arabs Theon's rule, \/a^ -\- b = a -^ , was used by 

2a 

Albategnius (c. 920), Abfll-Wefi (940-998), al-Nasavi (c. 1025), 

' Consult the introduction to C. J. Gerhardt's Das Rechenbuch des Maximus Planu- 
des, Halle, 1865. 

' Treutlein, "Das Rechnen in i6. Jahrhundert," in Abhandl. Math. Wissenschaft 
(1877); Triparty in Boncompagni's BuUeiino di bibliografia e di storia delle scienze 
malemaliche e fisiche (hereafter referred to as Boncompagni), XIII (1880), p. 695; 
Cantor, II (1900), p. 411. 

' Colebrooke, p. 145. 

* Kaye, Indian Mathematics, Calcutta, 1915, pp. 2-8; S. Gunther, Quadratische 
Irrationalileten der Allen, pp. 39-42; p. U2. 

' Colebrooke, p. 155, Example 51. 
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and al-Qalasidl (c. 1475).'" This rule was also used by the Jewish 
writer Johannes Hispalensis (c. 1140). Al-Karchl (c. 1020) used a 

b 



different formula, namely \/a^ + b = a -\ . How he de- 

2a + I 

rived this formula is not known ; possibly he used a method of double 

false position and interpolation by proportional parts similar to 

the one used later by Leonardo of Pisa for cube roots. Ibn AlbannS 

(c. 1300) used al-Karchi's formula for b > a, but Theon's formula 

for b ^ a. Leonardo used Theon's formula with a corrective sup- 



(-) 



plement, as follows: Va^ + b = a -\ — . A formula 

2a / , ^ \ 
2la + -) 

frequently used in the Middle Ages was V^ = - y/Aa"^. In the 

a 

Renaissance period Bombelli (1572), Cataldi (1613), and Schwenter 

(1618) approximated irrational square roots by the use of continued 

fractions. 

A formula for approximating irrational cube roots was invented 

by Leonardo of Pisa." In deriving this he used double false position 

and interpolation by proportional parts. He lets a' and (a + i)' 

be the two cubes nearest to a' + b. Then, since a'< a' + & < 

(a + i)', we have o^< 5 < 3a (a + i) + i. This relation gives a 

criterion for the proper choice of a. Also, if we increase a by unity 

we increase the resulting cube by 3a (a + i) + i. What increase in 

a will increase the resulting cube by 6? By proportion 

b 



he obtains as an answer . Therefore "i/d? -\- b = 

3a(a-fi) + 1 

a + . He gives these illustrations of his formula: 

30(0 + i) + I 

'" Gunther, Quadratiscke Irrationaliteten der Alien, pp. 45-46; Cantor. I (1908), 
pp. 633. 641; Hankel, Zur Geschichte der Malhematik im AUerthum und Mittelalter 
(hereafter referred to as Hankel), Leipzig, 1874, p. 185. 

" Boncompagni, Leonardo Pisano, I, pp. 370, 378, 380-81; Cantor, II (1900), pp. 
31-32; Bibl. Math. (3) II (1902), pp. 350-54. 

" Juan de Ortega (1512) used a modification of this formula, viz., •C'aa -|- j = 

b 
« + , , ,• See Cantor, II (1900), p. 388. 
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■^900 = 9 H = 9 - approximately, and v^2345 

271 3 



13 + 



148 

547 



= 13-, approximately. 
4 
To find the roots of large irrational numbers two plans were 
followed, the Hindu arrangement being used in both cases. Leo- 
nardo illustrates them both.^' In one of his solutions, V'927435 = 

\32I/ 



963 + 



II 



321 

follows : 



'-321 

2(962 + — ) 
\ 321/ 



the derivation of 963 is shown as 




We find this method used as late as the seventeenth century 
in the Kholasat-al-Hisdb " of Beh§. EddJn (c. 1600). Another plan 
is followed by Leonardo in this solution: 

V7234 = — V 72340000 = . 8505 - = 85 . 

100 1000 4 20 400 

This plan became increasingly popular and when in 1539 a man of 
Cardan's prestige adopted it systematically in his arithmetic,'* 
accompanied by clear and concise rules, it became the standard 
method in Europe. 

It had taken a thousand years, from Aryabhafa to Cardan, to 
perfect this approximation method. But it is the most effective 
and most widely used method known. The methods of Vieta and 
of Horner are its lineal descendants. 

1" Boncompagni, Leonardo Pisano, I, p. 355; Cantor, II (1900), p. 30; Curtze, 
"Ueber elne Algorismus Schrift des XII. Jahrhunderts" in Abhandl. Geschicht. Math., 
VIII (1898), pp. 1-28. 

" Taylor's Lilawali, Preface, p. 14; Hutton's Tracts II, London, 1812, p. 180. A 
manuscript copy of the Kholasat-al-HisSb is found in Professor D. E. Smith's private 
collection. 

" Cardan's Arithmetic, chap. 23. 
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3. SUMMARY 

1. The superior symbolism and notation of the Hindus enabled 
them to invent and use with ease formulas based on the principle 
of inversion. 

2. For finding the rational roots of numbers too large for hand- 
ling by inspection they used a method of exhaustion. It was based 
on the inversion of the expressions a^ + 2ab + b^ = (a + by and 
a" + 3a^6 + 306'' + 6' = (a + by, and on the decimal notation with 
its place values. Its essence was the evolving of the digits of the 
root by orders (Evolution). This method was also used for irra- 
tional roots later on, and European scholars standardized an adapta- 
tion of it for approximating irrational roots by decimal fractions. 

3. The five principal approximation formulas invented and used 
in the Middle Ages were: 



(o) Va^ + b = a+ -; (b) Va^ + b = a + 



2a 2a -+- I 



(c) VA = ~ \/Aa^; (d) ^a^ + b == a + 



a 3c(a -f i) -f I 



(e) ^a^ + b = a + 

30(0 + i) 



IV 

APPROXIMATIONS FOR SPECIFIC PURPOSES 

I. THE RISE OF NUMERICAL HIGHER EQUATIONS AMONG 
THE ARABS 

Although the Arabs did not contribute much original matter to 
algebra they vitalized it and enriched its contents by applying 
algebraic operations to the problems of Greek geometry and to 
their own problems in astronomy and trigonometry. This led them 
directly to numerical higher equations. Archimedes's problem of 
the section of the sphere, which led to a cubic, was first attempted by 
al-Mclh&nl (c. 860), and was later solved successfully by al-Kh&zin 
(c. 950) with the aid of conic sections. The trisection of angles and 
the computation of sides of regular polygons led to other cubics,' 
which were also solved by conic sections; for they soon came 
empirically to the conclusion that cubics could not be solved alge- 
braically. The solution of numerical cubic equations by intersecting 
conies was the greatest original contribution to algfebm made by 
the Arabs. These solutions remained unknown to the Western 
world, and were rediscovered in the seventeenth century by Des- 
cartes, Thomas Baker and Edmund Halley. The success of the 
Arab scholars in this field may have deterred them from trying 
methods of approximation. What they might have done in this 
field may be inferred rather than judged from the solitary example 
left us in Saracen writings. 

2. A SOLUTION BY AN UNKNOWN ARAB SCHOLAR 

In Miram Chelebi's annotations of Ulugh Beg's astronomical 
tables (1498) there is explained a method of solving the equations 
Px = x^ + Q, where P and Q are positive and x' is small in com- 
parison with Q.^ Chelebi explains that this equation was used for 
a specific purpose, namely, to find the sine of 1°. With the Persian 
astronomers, as with us, the ordinary trigonometric interpolation 

^ L'algibre d'Omar Alkayyami (translated and edited by Woepcke), Paris, 1851, 
pp. 54-57; p. 82; Hankel, p. 277. 

! Cantor I (1907), pp. 781-82; Hankel, pp. 289-93; Braunmiihl, Gesckichte der Vor- 
lesungen der Trigonomeirie I, Leipzig, 1900-3. PP. 72-74- 
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would fail for sin i°. The solution is commonly attributed to 
Giyit Eddin al-KSschl. Braunmiihl, however, thinks it is due to 
al-Zarkaii (c. 1080). 

In this particular equation P = 45' and Q = 47' 6* 8' 29", where 
the radius length i' is the sexagesimal unit and is equal to 60*. The 
solution is as follows: 
From the equation, 

+ x' Q 

X = ^-^ — ; for a relatively small x we have x = -. 

Then x = a + r;a + r = = a -] 

,R + a' 
= a+^. 

Suppose = b, with remainder S; that gives R = bP + S — a'. 

P 

Then x = a + b + s; a + b + s= ^^ 

,R + {a + b + sy ,R + {a + b + s)' 
= a + =a+ 

S-a'+ja + b + sy „ , , , S+ja + by- a' 

= a + b -\ = a + b -] 

P P 

Suppose the last fraction equals c, with a remainder T. 

Then x = a-\-b-\-c-{-s, and so on. 

The actual problem is to evaluate sin i" from the known trigo- 
nometric relation 60' sin 3° = 3.60' sin 1° — 4 sin' 1°. Letting 
s; = sin 1° and simplifying, the equation takes the form 45' x = 
x' + ^ • 60'' sin 3°. The Arabs knew the last term to equal 
47' 6* 8' 29" = Q. 

The following' is Chelebi's problem, carried out to seconds; in 
the Arab text it is carried out to fourths. 

' Hankel, p. 292. 
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Divisor P 


45' 


a 


b 


c 


Quotient x 




I* 


2' 


49" 


Dividend Q 


47" 


6^ 


8' 


29" 


aP 


45 








R 


2 


A' 




I 






R a' 




127 


bP 




90 






S 


37 


{a + by - a' 






6 


12 


(a + by - a? 




2234 


41 


cP 






2205 




T 


29 



That the method was effective for its particular purpose is shown 
by Ulugh Beg's tables. That it had its limitations is also obvious: 
it is laborious; and it is restricted to equations where P and Q are 
positive and x is relatively small. Yet the advance in the theory 
of equations was soon to remove the latter limitations; and refine- 
ments and modifications by men like Harriot, Raphson, and 
Halley might have reduced the former. Certainly this remarkable 
solution was the feat of some genius seemingly unsupported and 
isolated. 

We quote Hankel's appreciation : 

"This beautiful method of solving numerical equations stands 
second to none of the methods of approximation invented in the 
Occident since the time of Vieta in fineness and elegance. Apart 
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from the solutions of square and cube roots, with which it has a 
similarity of principle, it is the first method of successive numerical 
approximations which we meet in the history of mathematics." 

3. A SOLUTION BY LEONARDO OF PISA 

The last statement may or may not be true. There were at least 
two other modes of approximation in the thirteenth century. One 
was a well-developed system in China which we shall describe in 
the chapter on Horner's method. The other was the method em- 
ployed by Leonardo of Pisa. 

At a scientific tournament in Pisa in 1225 he was asked by the 
Emperor Frederick II to solve the equation x^ + 2x^ + io:>c = 20 
by Euclid's method of ruler and compasses.^ After several unsuc- 
cessful attempts he changed his method of attack and showed by 
rigorous analysis that the roots of this equation could not be found 
by ruler and compasses. Then he finally gave this approximate 
solution : X = 1° 22' 7" 42'" 33 '^ 4 ' 40 ". What process he used no 
one knows. But after seven centuries of further study in calcula- 
tion we differ from his result by onlv ij^" or 

31,104,000,000 

As with the solution of the Arab scholar, so with this one : we marvel 
alike at its isolation and at its accuracy. 

Many attempts have been made to explain Leonardo's method. 
Zeuthen thinks he used a method of false position with interpolation 
similar to his method for the cube root of numbers. Others think 
there was some relation between his work and that of the Arab 
scholar. Professor David Eugene Smith thinks he was acquainted 
with the work of the Chinese. There is evidence in his writings 
that he was familiar with the Chinese method of approximating 
roots of indeterminate equations.' 

4. SUMMARY 

I. The Arabs added little that was original to the science of 
algebra. But they vitalized it by applying the new science to 
problems of geometry, astronomy, and trigonometry. 

* Cantor, II (1900), p. 46; Boncompagni, Intorno ad alcune opere di L. Pisano, pp. 
17-19. 

' Cantor, II (1900), p. 26; Curtze in Zeitschrift fiir Malh. u. Phys. XLI, Hist. Lit. 
Abt., pp. 81, 82. 
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2. This gave rise to higher numerical equations, many of which 
they solved by intersecting curves. The success of this method 
undoubtedly delayed the work of solving them algebraically by 
approximation. 

3. We have only two extant algorithmic solutions by approxima- 
tion from this period, namely that described by Miram Chelebi 
and the solution of Leonardo of Pisa. They were both effected for 
a particular purpose. 



V 

METHODS INVOLVING DOUBLE FALSE POSITION 
AND PROCESSES OF INTERPOLATION 

I . THE GENERAL NATURE OF THESE METHODS 

The first successful general methods of approximating roots in 
complete equations were those involving "double false position," 
and the history of these methods is that of incomplete equations 
repeated. The chief problem in this mode of attack is to discover a 
practical way of finding intermediate ("mean") values that will con- 
verge to the true root. Three chief plans sprang up : the "averaging" 
rule of Chuquet; the methods of interpolation by proportional parts 
used by Cardan, Burgi, and Pitiscus; and the process of orderly 
evolution of digits used by Stevin, known later as Stevin's Rule. 

2. CHUQUET's RULE OF MEAN NUMBERS 

Chuquet's Triparty ' (1484) was little known outside of France. 
It was circulated only in manuscript, not being printed until 1880. 
Yet an examination of the textbooks of LaRoche, Petro Sanchez 
Ciruelo,^ and other writers reveals that he had considerable influence 
in France. His was a mathematical mind of a high order. His 
symbolism and arrangement of work were much ahead of his time. 
In the theory of equations he held advanced views; he recognized 
negative and hinted at imaginary roots, and gave a definite state- 
ment for the law of signs. 

In his Rule of Mean Numbers he states, without proof, that 

ffli + 02 ffli a2 

r — — - lies between — and rr. Also, if &i and &2 have the same 

Oi + O2 Oi O2 

signs, the following differences will have the same sign: 

a\ ai -j- 02 0162 + ffl2&i ai -\- Oi a-t aibi -\- ajbi 

bi bi + &2 bi{bi + bi) bi + bi b^ b^ (&i + 62) 

The method is perfectly general, and he applies it to complete and 

^Triparty, pp. 653-54 in Boncompagni, XIII (1880); Cantor, II (1900), pp. 351, 
359-60; Boncompagni, XIV (1881), pp. 416, 419. 
2 Cantor, II (1900), p. 386. 
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incomplete equations alike. An integer a is put in the fractional 

a 
form-. For rational roots the approximation terminates of itself; 

for irrational roots the rapidly increasing denominator insures a 
more rapid convergence than the arithmetic mean that had been 
used by others.' As an illustration we may take his solution of the 

equation »:■' -\- x = 39—. 

81 

In this case we see that :>(;=- is too small, while x = - is too large. 

I I 

Hence the first mean is ^ = — . By substitution this proves 

I + I 2 

too small. Therefore the new limits are — and -, the new mean 

2 I 

being = — . This second mean proves too small. The lim- 

2 + 1 3 

17 6 . . 17-1-6 23 
its — and - give the third mean — = — . This is also too 

31 3+14 

• 1 rr, ,. . 23 , 6 . 23 -1- 6 29 , , 

small. Ihe hmits — and - give = - — , which proves too 

4 I 4+1 5 

23 29 23 -1- 29 S2 

large. The limits — and — give the mean -^^ = — , and this 

4 5 4+5 9 

proves to be the exact root. 



3. THE REGULA AUREA OF CARDAN 

The first general method of approximating roots of numerical 
higher equations to be printed and publicly known was the "Regula 
Aurea" (Golden Rule) of Cardan (1501-1576) published in his 
Ars Magna ■• (1545), six years after the publication of his rules for 
square and cube roots. It is built on the basis of two false positions 
and a particular mode of interpolation. Though he uses it only for 
equations of the third and fourth degree, it is applicable to equa- 

' An adaptation, by no means an improvement, of Chuquet's method was worked 
out by DeLagny in his Method of Mediation (1697). 

' Chapter XXX, entitled "De regula aurea"; Cantor, II (1900), pp. 505-507. 
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tions of every degree. He restricts the rule to equations wliere the 
coefficients are positive or zero. 

The following is his derivation of the rule, presented in modern 
analytical form: 

Let/(x) = k, where /(x) is a polynomial in x arranged in descend- 
ing order. Let a and (a + i) be two positive integers such that 
/(a) = k — b and /(a + i) = k + b'. Then the root x lies between 
a and a + i. That is, x = a + t{a + i — a) or a; = (a + i) — 
e{a +1 —a), where the corrections t and e are positive proper 
fractions. The first may be called the additive, the second the sub- 
tractive form. Since /(a -|- i) > f(x) > f{a), we have f{a + i) — 
f{a) > /(x) — /(a) > o and consequently 

f{x)-f{a) _ k-{k-b) _ b ^^ 



f{a + I) - /(a) {k + b') - ik-b) b + b' 
Take as an approximate value of I. Then x = a- -\- 



b + b' b + b' 

Since < i our next additive form will then be/ ( a -| I 

b + b' \ b + b'J 

= k - b". 

Placing this new additive value in the subtractive form we have 

In the subtractive form e is less than unity. For an approximate 

value sete = f(-+^)-fM = jk + b') - k ^ 

/(a+x)-/(a+^-^,) (^+^)-(*-^") 

-. This gives the third approximation x = a + i — 



b' + b" ■ b + V 

b' 
b' + b 



- = f, say. According as k lies between f(r) and /(a) or be- 



tween /(r) and /(a -\- i) one must, in continuing, use the additive or 
the subtractive process. 
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His first illustrative example is the solution of the equation 

X* + T,x^ = 100, the root of which he finds to be 2 ^^. The fol- 

4697 

lowing is his solution and arrangement: 



2 

4io;- 



■100- 



60 



62 



122- 




162 



252 

61 



31 

61 



31 
61 

1922 

4697 



1922 

62 -^— 77 



61 



2775 
4697 



3 
-162 



Contrary to his own restrictions on the signs of the coefficients 
Cardan uses the rule in the equation x^ = 6x + 20, with disastrous 
results. 



4. stevin's contribution 

The approximation work in Stevin's algebra (1685) might with 
equal propriety be discussed in connection with Vieta or Newton. 
For Stevin (1548-1620) anticipated Newton's analytical form and 
treatment as well as Vieta's idea of evolving the digits of the root 
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of a complete equation order by order. But since his scheme of 
deriving the digits involved the successive choice of two positions, 
we discuss his method in this chapter. We can best explain his 
process by recording one of his illustrations.^ 

Suppose x' = 300X + 33915024; 
ii X = I, then i^ < 300- 1 + 33915024; 
if:*:= 10, then io3< 300-10 +33925024; 
i{ X = 100, then 100' < 300=100 +33915024; 
ii X = 1000, then 1000' > 300-1000 + 33915024. 

Hence x lies between 100 and 1000. By trying x = 200, x = 300, 
X = 400, he finds the root to lie between 300 and 400. Similarly 
trying 310, then 320, then 330, he finds x to lie between 320 and 
330. A similar procedure for the units gives 324 as the exact value 
of the root. Stevin observes that an irrational root can be approxi- 
mated to within any desired degree of accuracy by using this 
method and his scheme of computing decimal fractions. This 
method was advocated and used considerably by Albert Girard in 
his Invention Nouvelle en I'Algebre (1629). Though laborious, the 
method is general. "Stevin's rule" was used by later algebraists, 
like Oughtred (1631), Kersey (1673), and Saunderson, in connec- 
tion with other methods.^ 

5. THE METHODS OF BURGI AND PITISCUS 

Pacioli's Summa (1494), the earliest printed work on arithmetic 
and algebra, gave to Renaissance scholars the method of double 
false position and interpolation by proportional parts used by 
Leonardo for finding roots of numbers. Leonardo's method now 
came to be employed frequently also in the solution of higher 
numerical equations, and Pacioli even used it in an exponential 
equation.' This method, which has always been basic in finding 
trigonometric functions, was used by Jost Biirgi in solving equa- 
tions which rose out of his work in trigonometry.' We give his solu- 

' Stevin, Les CEuvres Mathematiques (editor Girard, 1634), I, problem 77; Cantor, 
II (1900), p. 628. 

' N. Saunderson, Elements of Algebra, II, Cambridge, 1790, p. 728. 

' Pacioli, Summa, fol. 186 recto to 188 recto; Cantor, II (1900), pp. 325-26. 

' Cantor, II (1900), pp. 626, 645-46, 648; C. I. Gerhardt, Geschichte der Math, in 
Deutschland, Munich, 1877, pp. 76-82; 84-87. 
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tion of such an equation, namely, 9 - 30*;^ + 27a;* — ^x* -\- x^ = o. 

The root has been found to lie between 0.68 and 0.69; let x = 0.68 

+ h, and f{x) = o. Then 

/(0.68) = 0.0569 /(0.68) = 0.0569 
/(0.69) = - 0.0828 /(0.68 + h) = 0. 

Difference o.oi : 0.1397 Difference h:o .0569 

Hence 0.1397 '■ 0.0569 = o.oi -.h; or h = 0.004, and x = 0.684. 
Repeating the process with /(0.6840) and /(0.6841) he gets 
0.00140012 : 0.00056410 = o.oooi : h; or h = 0.00004029. 
Hence x = 0.68414029. 

A more conventional method of double false position is used by 
Pitiscus (1612), a pupil of Biirgi. In solving the equation' 

^3 
5,176,380 = 30: — a;' he sets x = 1,725,460 H . Manifestly x > 

3 

1,725,460. By trial he finds x to lie between 1,730,000 and 1,740,000. 

Let us call these two values ai and a^ and represent the equation by 

fix) = k. Then/(ai) = k — 38,158 = k — dufia^) = k — 9,061 = 

k — di. Then 

Oidi — aidi 

X = — -— = 1,743.1 14- 

d\ — di 

6. SUMMARY 

1. The first successful general method of solving complete 
numerical equations by approximation was given by Chuquet in 
his Triparty (1484). It was not printed until 1880 and remained 
comparatively little known. 

2. The first printed method was Cardan's Regula Aurea (1545). 
It was restricted as to the signs of the coefificients but not as to the 
degree of the equation. It never came into general use. 

3. In 1585 Stevin explained a method of evolving the digits of 
the roots by order. 

4 . Pacioli and Biirgi extended Leonardo's method of finding cube 
roots to the solution of complete equations of higher degree. This 
method and the conventional method of double false position ex- 
emplified by Pitiscus were the ones most generally used before the 
year 1600. 

• Pitiscus, Trigonometriae, pp. So-53. 
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VIETA'S EXTENSION OF THE HINDU METHOD 
OF APPROXIMATION TO COMPLETE EQUATIONS 

I. VIETA's METHOD OF APPROXIMATION 

A new epoch in the history of the solution of numerical equations 
begins with the year 1600, the date of the first publication cf 
Vieta's De Numerosa Potestatum Purarum atque Adfectarum ad Exe- 
gesin Resolutione Tractatus. In it was explained a systematic and 
concise way of solving any equation, complete or incomplete, by 
methodically finding the successive digits of the root, beginning 
with the highest order. Why no one before Vieta should have 
thought of applying to the solution of complete equations the 
Hindu method of finding roots of large numbers may seem strange 
to us unless we reflect that useful inventions usually seem very 
natural after they are invented. At any rate, this method had been 
applied to large numbers for a thousand years, had appeared in 
definite forms for four hundred years, and had been perfected in its 
finest details, before any one thought of extending it to complete 
equations. It was the first comprehensive method of solving such 
equations that had been attempted, and it involved no restrictions 
as to terms, signs, or degree, for it applied to a form into which 
any equation can be changed. Vieta calls an equation "duly pre- 
pared" if it has the form x" -\- ax'^ + bx"'^ -\- = K where 

K is positive and all the coefficients are integers. 

De Numerosa Potestatum was first printed in Paris in 1600 
for private circulation. In 1646 Franz van Schooten, of Leyden, 
again printed it in a one-volume collection ^ of Vieta's works. 
It contains sixty-six pages, of which the first ten deal with in- 
complete equations. He arranges the solution in a way slightly 
different from the conventional one in order the better to show 
the analogy between the pure and the affected powers.^ We 

' Vietae Opera Malhematica, Leyden, 1646; hereafter referred to as Vieta. 

' Vieta classified powers and equations as pure and affected ("ad-f ected") . The 
names are geometrical in origin ; for Vieta considered all the terms in an equation homo- 
dimensional. Thus in the "affected" equation *' + 30* = K, the "resolvend" K is 
considered three-dimensional; and 30 is the "plane coefficient" in the "adjoined solid" 
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insert the solution for one pure equation and two types of afifected 
equations,' viz. : 

(i) x^ = 157,464; where a; = 54. 

(2) x^ + 30a: = 14,356,197; where x = 243. 

(3) 95,400a; -\- x^ = 1,819,459; where x = 19. 

The arrangement is that of Vieta. The explanations, too, are a 
close paraphrase in modern notation of his marginal and other 
explanations rhetorically expressed. 

(i) Let /(a;) = K for x> = 157.464; 

lix = a -\-b, then a' + sa^ft + 306^ -\- b^ =. 157,464. 
"Proponantur i C, aequari 157,464. 

I. Eductio lateris singularis primi 



K : Resolvend 

a' 

K -fia) : Residual 



157 


464 


125 




32 


464 



(a 



32 


464 


7 


5 




15 


7 


65 


30 




2 


40 




64 


32 


464 



II. Eductio lateris singularis secundi 
K -f(a) : Resolvend 
3a8 
30 
Divisor 7 65 (b 

sab' 

Subtrahend 

Itaque i C, aequetur 157,464 fit i iV 54." 

(2) Let /(a;) = K for x' + 30a; = 14,356,197; 

\l x = a + b, then (a' + za^b + saft^ + b^) + 3o(a + fe) = 
14,356,197- 
"Proponantur i C + 30 iV, aequari 14,356,197. 

30a:. The names "pure" and "affected" continued in use until the latter part of the igth 
century. Since, in the modem theory of equations, the basic form is aox" + aix" ~ ' + 
... + K = o, recent writers use complete and incom-plete as being more descriptive. 
In discussing the works of Vieta and his immediate successors, we shall employ the 
terminology used by them. 

'Vieta, pp. 167-68; 177-78; 179. 
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8 


6 





8 


oo6 





1 6 


350 


197 
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I. Eductio lateris singularis primi 

Plane coefficient 

K: Resolvend 

a' 
300 

m 

K-f{a): Residual 

II. Eductio lateris singularis secundi 

Plane coefficient 
K -f{a) : Resolvend 

30' 

3a 

Divisor 

3a"& 
3a&2 

h' 

30b 

Subtrahend 

K-f{ab) 



(a = 2 







3?. 


i 6 


350 


197 


I 


2 
6 




I 


260 


30 


4 


8 

96 
64 






I 


20 


5 


825 


20 




524 


997 



(& = 4 



III. Eductio lateris singularis tertii 





30 


524 


997 


172 


8 




72 


173 


550 


518 


4 


6 


48 




27 




90 


524 


997 



(6 



Plane coefficient 
Resolvend 

3a 

Divisor 

Sa^b 

Sab^ 

b^ 

306 

Subtrahend 

Itaque i C30 iV, aequetur 14,356,197 fit i iV 243." 
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(3) Let/(a:) = K for 9S,400x; + x;' = 1,819,459; 

\l X = a -\-h, then 95 400(0 + b) + {a^ + 30*6 + 306* + 6') 
= I-8I9-459- 

"Paradigma cum solidum adfectionis sub latere majus 
est cubo. 



9 


540 





I 


819 


459 





954 


00 


I 


819 


459 




954 

I 


00 




955 






864 


459 



I. Eductio lateris primi inanis ante devolutionem 
Plane coefficient 
K: Resolvend 

Quoniam 9 major est unitate, fit devolutio. 

II. Eductio lateris singularis primi post devolutionem 

Divisor 954 00 (a = i 

K: Resolvend 

95 400a 
a' 

/(a) 

K-f{a): Residual 

III. Eductio lateris singularis secundi 
Plane coefficient 
K -f{a) : Resolvend 
3a« 

31 

Divisor 95 730 (b = 9 

95 400J 

3a'b 

Sab^ 

b' 
Subtrahend 

Itaque si 95,400 N + iC, aequentur 1,819,459 fit 
I N 19." 

In a similar fashion he illustrates the solution of equations of the 
fourth, fifth, and sixth degrees. Vieta's work involves mainly four 
points: 



95 


400 


864 


489 




3 




3 


95 


730 


858 


600 


2 


7 


2 


43 




729 


864 


459 
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1. The first step is to ascertain the number of digits in the root. 
For pure powers, of course, this is easy. For affected powers there 
will also be as many digits in the root as there are periods in the 
resolvend if the highest power is dominant; that is, if the pure 
power is larger than the adjoined solid or solids. This will generally 
be indicated by the size of the coefficients. Thus examples (i) and 
(2) will have roots of 3 digits. A more difficult situation arises if 
the highest power is not dominant. Thus in example (3) the re- 
solvend has 3 periods, yet the root has only 2 digits. Such a case 
he calls "devolution." If the cube and the adjoined solid had been 
of opposite signs there would have been more digits in the root than 
periods in the resolvend. Such cases are called "anticipations." 
Vieta and all his followers found these cases full of practical diffi- 
culties and he offers many suggestions for reducing and circum- 
venting them. One device is shown in our third illustration. 

2. With Vieta the fundamental operation for finding the roots 
of pure as well as affected equations is the expansion of the powers 
of a binomial.* 

3. For each approximation a, Vieta takes the residual K -f{a) as 
the new resolvend. The entire residual is put down each time. 

4. Vieta gets his trial divisors by taking the sum of all the co- 
efficients of the powers of b except the coefficient of &", which is 
invariably unity. Both Cantor and Hankel have given erroneous 
explanations of Vieta's divisors. The mistake seems to come from 
representing Vieta's process in the compact equational form of 
Newton's scheme. At least they attribute to him Newton's divisor. 
Following Hankel's and Cantor's explanation, the first divisor in 
Example (i) should be 120,030. Vieta's divisor, however, is 
120,630.* In this connection it should be said that both these his- 
torians give an erroneous impression of the plan of Vieta's method. 
They lead one to suppose that Vieta's method of evolving digits 
is the same as Newton's method of substituting values in equa- 
tions of series; and Hankel says as much. But the two differ in 
origin and form, and in certain of the steps. 

In developing his approximation method Vieta uses for con- 

* Cantor's exposition of Vieta's process, II (1900), pp. 640-41, where he sets x = 
a:i + a;2 -t- X3, is apt to be misleading. Vieta did not use that substitution, though his 
expositors Hume (1636) and Herigone (1642) used it occasionally. 

' Hankel, p. 370; Cantor, II (1900), p. 641. 
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venience equations having rational roots. But he adds« that irra- 
tional roots can be approximated in the same way after first trans- 
forming the equation into one whose roots are decimal multiples 
of the roots of the given equation. This, of course, is nothing but 
an extension of Cardan's rule for irrational numbers. Thus, in the 
equation x' = 2, the root is found from solving x' = 2,000,000, 
; likewise if x' -|- 6oox = 8,000, the root correct to hun- 
dredths is found by solving x' + 6o,ooox = 8,000,000. 

2. THE MODIFICATIONS OF VIETA's METHOD BY HARRIOT 

It was not in his homeland that Vieta's method met with the 
greatest appreciation. The people of the British Isles have always 
been partial to mathematical processes dealing with numerical 
calculation. Three English scholars, Thomas Harriot (1560-1621), 
William Oughtred (1574-1660), and John Wallis (1616-1703), all 
noted as teachers and writers, became the chief protagonists of the 
new method. 

In his Artis Analyticae Praxis ' (written before 1621 and pub- 
lished posthumously in 1631) Thomas Harriot introduces two devices 
for expediting the calculation in Vieta's method. One of these is 
the elimination of the second term; the second is his Canones 
Directorii, a collection of the different cases of numerical solutions, 
ready and arranged for use under various canonical equations, 
showing the elements necessary to form the several subtrahends and 
resolvends. This plan gained much vogue and we find Raphson, 
Halley, and DeLagny doing the same thing in connection with 
Newton's method. 

In the actual computation he uses Vieta's method unchanged. 
We find no basis for De Morgan's statement that he formed "only 
so much of the divisor as is necessary for the determination of the 
next figure." Nor does WaUis, who is very partial to Harriot, men- 
tion that fact. We subjoin the first two steps of an example from 
page 156 of the Praxis, which is quoted in full by De Morgan. 

' Vieta, p. 228. 

' The writer has been unable to consult the original and depends for his description 
on the following secondary sources: De Morgan's article, "Notices on the Progress of 
the Problem of Evolution," in Companion to the British Almanac, 1830, pp. 34-52; 
Hutton's Tracts, II, pp. 278-86; Cantor, II (1900) ; Henry Stevens, Thomas Hariot and 
his Associates, London, 1900; J. Wallis, A Treatise of Algebra both Historical and Prac- 
tical, London, 1685. 
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"Sit aequatio numerose proposita 

aaaa - 102400 + 6254a = 19633735875 

Radix universalis successive educenda. 



Homogeneum resolvendum 



hhhh 



3 7 5 
19633735875 







ggg 


6254 






-If 


— 1024 






bbb 


-27 






+ 


27006254 






Divisor A 


26903854 






gggb 


18762 


Radix sing, prima 6 = 3, 




-ffbb 


— 9216 






bbbb 


81 






+ 


81018762 




Ablatitium Ab 


80097162 




Radij 


: singularis 


3 


Homogeneum residuum resolvendi 


um 


11624019675 






ggg 


6254 






-If 


— 1024 






-2ffb 


— 61440 


Radix sing, decupl. b = 30, 




4bbb 


108000 






6bb 


5400 






46 


120 






+ 


1 13526254 






— 


62464 






Divisor B 


112901614 






gggc 


43778 






-ffcc 


- 50176 






-2ffbc 


— 430080 


Radix sing, secund. c = 7 




^bbcc 


756000 






6bbcc 


264600 






^bccc 


41 160 






cccc 
+ 


2401 




1064204778 






— 


480256 



Ablatitium Be 1059402218" 
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Harriot's canonical form is aaaa — ffaa + ggga = hhhh. He lets 
a = b + c. Then the transformed equation becomes {bbbb + 
4.bbbc + 6bbcc + ^bccc + cccc) — // {bb + 2&c + cc) + ggg {b + c) 
= hhhh. By Vieta's method the divisor would be {ggg — // — 
2//6 + 4666 + 6&6 + 4i>) ; and this is identically the divisor in the 
illustration quoted, both in the second and the third step. 

Harriot's Praxis was the most analytical algebra that had been 
written up to that time. His lucid treatment of the theory of equa- 
tions, his improved symbolism, and the inspiration of his person- 
ality to his students were accountable in no small measure for the 
popularity of Vieta's method in England, and it was often called 
after him "Harriot's method." 



3. THE IMPROVEMENTS OF OUGHTRED AND WALLIS 

No one did more to popularize the new method than did the 
clergyman mathematician William Oughtred. This he accomplished 
by giving private tuition to ambitious young men and these spread 
his teachings throughout Great Britain; among them were Seth 
Ward, Christopher Wren, and John Wallis. An even more potent 
influence was his textbook, the Clavis Mathematica. Before the 
year 1700 five editions of this little volume had gone through the 
press. 

The first edition appeared in 1631, the same year in which 
Harriot's Praxis was published. The Clavis appeared first. In it 
there is no reference to Harriot, nor any use made of his improve- 
ments : the two men evidently knew nothing of one another's work. 
Oughtred uses Vieta's method throughout. In "locating" the root 
he depresses the equation and uses a method similar to that of 
Stevin for obtaining the first digits of the root. It is interesting 
to observe that he uses the decimal bar and logarithms in his 
solutions. 

John Wallis made this method more practical by arranging the 
work in a compact, yet clear form, so as to show the steps in the 
computation. He differed from Vieta, Harriot, and Oughtred in 
not rejecting the coefficient unity in his divisor. To help discover 
cases of devolution and anticipation he employs two sets of periods, 
one for the resolvend and another for the coefficient; at times he 
uses Stevin's rule for this purpose. 
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4. THE CONTRIBUTION OF HUME 

In France the method of Vieta made slower progress. In 1636, 
James Hume, an emigrated Scotchman, published in Paris a work 
with the title Algebre de Viete, d'une methode nouvelle, claire, et facile. 
It is an abridged presentation of Vieta's analytical works, repro- 
ducing both his methods and his illustrations. However, he intro- 
duces one modification : he uses only so much of the divisor and of 
the dividend as is needed to compute the next digit. Thus in the 
example (2) which we quoted from Vieta he uses 12 for the divisor 
and 6350 for the dividend, where Vieta uses 126,030 and 6,350,197.* 

Vieta's method was adopted by P. Herigone^ in his textbook 
Cursus Mathematicus, published about this time ; but neither Hume's 
nor Harriot's improvements were used. 

5. SUMMARY 

1 . The Hindu method of finding the roots of large numbers was 
extended by Vieta to the solution of complete equations in his De 
Numerosa Potestatum (1600). 

2. The method was first taken up in England by Harriot (some 
time before 162 1), who modified it by using canonical equations in 
which the second term had been eliminated. 

3. Oughtred (1631) and Wallis (1685) improved the arrange- 
ment, but otherwise used Vieta's method unmodified, except that 
Wallis included the coefficient unity in his divisor. 

4. In France James Hume (1636) used only so much of the 
dividend and divisor as was necessary to evolve the next digit of 
the root. H6rigone's algebra (1642) was the first French text- 
book to popularize Vieta's method in France. 

' Hume, Algebre de Viete, pp. 234-38. 

' Herigone, Cursus Mathematicus, II, Paris, 1642, pp. 270-76. 



VII 
NEWTON'S METHOD OF APPROXIMATION 

I . ADVANCE IN THE THEORY OF EQUATIONS 

The seventeenth century saw a remarkable advance in the 
theory of equations. The symbolic notation introduced by Vieta and 
improved by succeeding algebraists, notably by Harriot and Des- 
cartes, made the symbolic equation a chief instrument of analysis. 

When Harriot put the equation in the form ao*;" ' + ai^c""' + 

+ K = o, the eye began to aid the mind in grasping the relations 
of the parts of an equation. The law of signs, used by Harriot, 
enunciated by Descartes, and refined by Newton, gave an effective, 
though limited, criterion for the number and nature of the roots 
in an equation. The invention of coordinate geometry by Descartes 
in 1637 and the independent discoveries of the differential calculus 
by Newton and Leibniz a few years later gave a powerful impetus 
to work in theory of equations and in practical numerical calcu- 
lation. 

For purposes of numerical approximation one of the tools made 
most effective by the improved symbolism was the infinite series 
in an equation. Series, as a means of exhaustion and of summation, 
had been successfully used for approximating numerical values 
since the time of Archimedes. But with the introduction of sym- 
bolic algebra, series gained an effectiveness undreamed of by the 
ancients. Pioneers in the use of equations in series for purposes of 
numerical approximation were Wallis (1656), Lord Brouncker 
(1620-1684), Gregory (1667), and Newton (1669). 

2. A DESCRIPTION OF NEWTON'S PROCESS 

It was while finding the area under curves by his newly-invented 
"method of fluxions" that Isaac Newton (i 642-1 727) discovered a 
new and most necessary application of series. Following the lead 
of Wallis he transformed such fractional and radical expressions as 



y = and y = Va^ + x^ into integrable form by expanding 

b + X 
them into equations in series by ordinary division and finding of 



34 Approximating Roots of Numerical Higher Equations 

roots. This scheme he extended to affected equations. "Since the 
whole difficulty lies in the solution," he says, in explaining it, "I 
shall first illustrate the method I use, in a numerical equation." ' 

Such then is the historical setting of the beginnings of a type of 
approximation by successive substitutions that is almost universally 
known as "Newton's method of approximation." It was entirely 
incidental to his work in quadrature by calculus. The originator of 
this famous method gives only a single numerical example, and it 
is illustrative and introductory to the more general problem of 
deriving a method for finding the root of literal equations by means 
of equations in series. Some historians notwithstanding, it is 
hardly correct to say that it was a further development or an adap- 
tation of Vieta's method. Of course, Newton's method includes 
Vieta's ideas of limit, substitution, and exhaustion; but these seem 
inherent to any method not using double false position. In some 
of its features Newton's solution resembles the approximation 
processes in Chelebt's annotations and of Stevin. Later algebraists^ 
sometimes used "Stevin's rule" of substituting i,io,ioo etc., to find 
the integral part of the root. Newton himself says of this method : 
"Whether this method of solving equations be commonly practised 
or not, I cannot tell; but surely to me it appears simple in com- 
parison with other methods and more suited to practise."^ 

Newton's solitary numerical example, the solution of y' — 2y — 5 
= o, has become classic. Lagrange, RjUffini, and Horner, all demon- 
strated the virtue of their methods by applying it to Newton's 
equation, and it has become a favorite problem in text-books. 

This process of approximating by successive substitution in equa- 
tions is thus explained by Newton: 

Suppose that an approximate value y = 2 has already been 
found, dilTering by not more than one tenth of itself from the true 
root. Then y = 2 + p. Substituting this value in the original 
equation gives />' + 6p^ + lop —1=0. Since p is a. fraction we 
can reject all but the linear terms, or lop —1=0 approximately. 
That is /) = 0.1, approximately; or ^ = o.i + g. Continuing this 
process as far as desired we get several values q, r, s, etc. And 

' Newton, Analysis by Equations of an infinite Number of Terms, translated and 
edited with commentaries, by John Stewart, London, 1745, pp. 321-28, 357. 
' N. Saunderson, Elements of Algebra, II, Cambridge, 1790, p. 728. 
' Newton, Analysis by Equations, p. 330. 
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y = p-^g.-\-r-\-s-{- If at any time the coefficients are 

such that one doubts the accuracy given by the Hnear terms, one 
may advantageously use the quadratic part of the equation. Thus 
6p'^ -\- lOp — I = o would give a more accurate value of p than 
lOp —1=0. In such cases, use the lesser root of the quadratic. 
Newton sets up his solution in paradigm form, as follows:* 



3,3 _ 


-2y-5=o 


+2,10000000 
-0,00544853 
+2,09455147 =3» 


2+p = y 


+3,3 


+ S+l2p+6p'+p^ 




-\-2y 


-4 -2p 




-5 


-5 




sum 


-i + iop+6p^+p' 


0, l+q = p 


+P' 


+0,001 +o,03g+o,3g2+g» 




+6p^ 


+0,06 +1,2 +6,0 




+ iop 


+ 1, +10 




— I 


-I, 


sum 


+0,061 + 1 1 ,23g+6,3g2+g3 


-o,oo54+r = g 


+6,3S^ 


+0,000183708 — o,o68o4r+6,3r2 




+ 11,232 


— 0,060642 +11,23 




+0,061 


+0,061 




sum 


+0,000541708 + 11,16196^+6,3?-^ 


— 0,00004854+5 = r 







He next takes up the more general method of literal equations.' 

"These things being thus shown in numbers, let y^ + a^y — 2(i' + 

axy — re' = o be proposed to be solved." Using the same method, 

X x^ 

settings = fl + />,hederives p= + 2; 2 = — +''> etc., and 

4 64a 

y = a--+ — + l^ + -5^ + A practical difficulty 

4 64a 5i2a^ 1638a' 



* Newton, Analysis by Equations, pp. 328-329. 
' Newton, Analysis by Equations, pp. 330-343. 
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in getting the supplementary approximation is encountered when 
the coefficients are letters. He reduces this by arranging the terms 
according to their dimensionality; this arrangement has been called 
Newton's Parallelogram. The value of y will converge when coeffi- 
cients are so chosen that x is small; and the smaller they are, the 
more rapidly the series converges. Later algebraists sometimes 
found numerical approximations by using as forms such solutions of 
literal equations. But Newton, as might be expected, did not sub- 
stitute numerical values of x until after the integration. With him 
the objective was the area of the curve and not the numerical 
value of y. 

Newton's description of his method was contained in an unpub- 
lished tract entitled De Analysi per aequationis numero terminorum 
infinitas which was shown to Barrow and Collins in 1669.^ His 
method was first given to the public in the 1685 and 1693 editions 
of Wallis's algebra. In 1704 and 171 1 it was published in Newton's 
own works in Latin and in 1736 and 1745 it appeared in Colson's 
and Stewart's translations. Newton's name was always linked with 
his particular procedure, but the appelation "Newton's method" 
came to be applied to this whole type of approximations through 
the writings of Lagrange and Fourier. 

3. raphson's development of Newton's method 

The first mathematician to enlarge upon the new method of 
approximating roots by substitution in equations of infinite series 
was Joseph Raphson. We may even say that he was the first to 
expand it into a system. It was he who made Newton's method 
a component part of modern algebra; for with Newton it was an 
accessory to the calculus, and was never taken up in his Universal 
Arithmetic. In 1690 Raphson published a tract of forty-eight pages 
entitled Analysis Aequationum Universalis seu Ad Aequationes 
Algebraicas Resolvendas Methodus Generalis, & Expedita, Ex nova 
Infinitarum Serierum Methodo. In 1697 he published a second edi- 
tion with an appendix discussing the recently published methods 
of DeLagny and Halley. Though quoting from Wallis ' as to 

' Letter to Mr. Oldenburg, Oct. 24, 1676, given in Stewart's commentary on Analy- 
sis by Equations, p. 358. 

' Raphson, Analysis Aequationum Universalis, London, 1697, p. 50, 
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Newton's work in infinite series, he does not even refer to Newton's 
method of approximation. 

Raphson's method is eclectic. He sets up the example like Vieta, 
employs Harriot's notation, and divides the coefficients and resol- 
vends into periods, "modo, quem ducere Vieta, nostratisque Har- 
riotus & Oughtredus." Like Oughtred he suggests the use of loga- 
rithms in finding the first digits of the root ; he adds that it may also 
be done by constructing parabolas and by examining the limits of 
the roots. But the main contributing elements in his method are 
Harriot's use of canonical forms and Newton's method of substi- 
tution. Employing and enlarging upon these ideas Raphson builds 
a canonical form for the supplement of the root in each of sixty- 
four canonical equations; so that instead of substituting anew 
in each equation, as in Newton's example, he could substitute the 
successive approximations in one canonical form. The "Newtonian 
divisor" in these forms is a misnomer and should be credited to 
Raphson. The name given to the collection of forms, "Canones 
Directorii," is the same as that used by Harriot. 

Only the two equations a^ = d and ba — a^ = c, representing 
pure and affected powers respectively, have the canonical forms 
for their supplements derived in full. Their derivations are as 
follows: ' 

(i) If a' = d, and a = g + x, then we have g^ + sg^x + 
:igx^ + x^ = d. Rejecting the powers of x above the first, as being 

d — g^ 
comparatively small, we have the supplement x = — . 

(2) If ba — a^ = c, and a = g + x, then "e Theoremata Vietae," 
i,g _ g3 ^ (J _ 2ig^)x — ^gx^ — x' = c. Again, rejecting the powers 

of X above the first, we have for the supplement x = — ■ — * 

b - U^ 

Then he proceeds to prove the convergency of the series. 

We append the following problem « solved by this last form : 



"Numeris 77284 a — aaa = 801 312 8 

' Raphson, pp. 5-6, 8. 
' Raphson, p. 19. 
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Theor. X = '^ + ^^^ " ^' 
b - 3Sg 

g = 200 

c + ggg-bg= +626328 

+ b - 3gg= - 42716) + 626328(- 14 = X 

200 

-14 

g = 186 

c + ggg-bg = + 143 160 

+ b - 3Sg = - 26504) + 143 i6o(- 54 = X 

186 

-5-4 

g = 180.6 

c + ggg-bg= +16152.216 

+ b - 3gg= - 20565.08) + i6i52.2i6(- 7854 = X 

180.6 

-7845 

g — 179.8146 vel 179.814 = g 

c + ggg - bg = + 322.299 405 144 

+ b - 3gg = - 19715.224) + 322.299405i44(- .016348 = x 

179.814 

— .016 348 



a = 179.797652." 

The tract contains the solutions of thirty-four numerical problems 
taken mostly from geometry and trigonometry. One problem, in- 
deed, deals with interest and is undoubtedly the first instance of 
Newton's method being used in the solution of exponential equa- 

J _ „365 

tions.'" Using* = — as the canon form for the supplement 

3652'" 
of the equation a'°* = b, he solves quite accurately the equation 
a^'^ = $1.06. 

4. THE APPROXIMATIONS OF DELAGNY 

Thomas Fantet DeLagny (1660-1734) spent over forty years in 
work on the approximation of roots. Observing that the expression 

'" Raphson, p. 36, problem 29. 
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z = y/a^ + b has a rational approximation z = a -\- — and a radi- 



2a 



■? f A 



cal approximation 3 = - + /f/ 1-&, he sought to obtain similar 

expressions for z = -v^a^ + &. These he found tohe - -\- A/ — h — > 

'? V A 



4 30, 



5' 'A 



b — I , ab + a ^ ^ , ab 
+ and a H > 2 > a + 



4 3a 3a2 + 6 + I Sa'' + & 

These results were published in 1691 '^ and in 1697 he incorporated 
them in the third part of his Nouveaux Elemens d'Arithmetique ei 
d'Algibre, and later in his Analyse Generale (1733). He derived 
similar expressions for higher powers. The process itself we shall 
discuss in connection with Halley's general formula. But the 
reader will be interested to know that in finding the roots of large 
numbers, DeLagny preferred to use the formula -s/a^ + b = 

a -\ rather than the Hindu method. Thus he finds the cube 

3ffl' + b 

root of 817,400,375 as follows: •v^8i7 = •v^729 + 88 = 9 ^-; 

multiply this by 100, since the desired root has three digits; the 
result is 935. 

DeLagny's method of attack in solving both pure and affected 
equations may be seen from his solution of the equation x' = aH + &." 
Let X = a -\- y. Then x' = a' + 30^ + 2>ay^ + y^\ also x^ = 
a» + a^y + b. Therefore b = ■^c^y'^ + 20^;)' + y'- Neglecting y' as 



being very small, 303'^ + 2a^y — b = o, or y = 1"4/ — ' ' 

3 '^ 9 3a 



Consequently, x = — + A/ 1 . Replacing & by & - y' in the 

3 ^ 9 30. 

" The incorrect data given in Cantor, III (1898), p. 115, seem to be taken from 
Halley's article in the Philosophical Transactions for 1694. 

" DeLagny, Nouveaux Clemens, Paris, 1697, pp. 303-6- Analyse Generale, Paris. 
1733. PP- 336. 342-46, 362-94. Journal des Scavans, May 14, 1691, pp. 297-98. 

" DeLagny, Nouveaux ^llmens, pp. 473-74- 
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last expression gives x = h 4/ 1 • ^^^ formulas for x 

3 '^ 9 3a 

and y are canonical expressions for the successive roots and supple- 
ments of the equation x' = a^x + h. Raphson's canonical form in 

c + 6? -f- ?' 
the equivalent equation a^ — ha = c is a = g + —^ 

3g^ - h 

different in form and derived by a different process. In DeLagny's 
illustration, x^ = 7569X + 240,903; then a^ — 7569, a = 87, & = 
240,903; this gives x = 58 + V1764 = 100 for a first approxima- 
tion, X = 58 4- /|/ 1755 for a second approximation, and so on. 

' 261 

DeLagny's "Method of Mediation" resembled the more effective 
Rule of Mean Numbers given by Chuquet; except that DeLagny 
used the arithmetical mean." 

Two points about DeLagny's work are worthy of special men- 
tion. He stood alone among his contemporaries in his insistence on 
a system of gauging the error of the approximation. "Every method 
is useless," he says, "if it be not accompanied by another method 
which gives the limits of error, be it par excis or par defaut." " And 
it was his pride that in the Triangle du Rapport he had an unfailing 
scheme for detecting the limits of error. 

The other point relates to his ambitious scheme of facilitating 
the solution of higher degree equations by means of tables. For this 
purpose he worked out three elaborate sets: (i) tables for solving 
nimierical equations directly, covering forty-nine large pages ; (2) the 
"Trapese Logarithmique," " for facilitating by means of arithmetical 
progression the solution of equations; (3) the "Triangle du Rapport" 
for expanding the roots into series. These tables together with 
Raphson's Canones Directorii might have come to fill a practical 
need, as do our integral and logarithm tables of to-day, had they 
not been made obsolete by the discoveries of Halley and Taylor. 

DeLagny's style was burdened with much verbiage and his 
classifications and arrangement lacked sharp lines. That may 

" Nouveaux Elemens, pp. 507-11. 

" Analyse GineraU, p. 546. 

" Analyse Generals, pp. 320-323; see also Memoires de VAcadimie for 1705 and 1706. 
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account somewhat for the comparatively slight influence of his 
voluminous work, some parts of which are very scholarly. Also, if 
he had lived in England instead of in France, his influence might 
have been more widespread. At any rate, the forty-eight pages of 
Raphson left a greater immediate impress on approximation 
methods than did the twelve hundred pages of DeLagny. 

5. halley's invention of general formulas 
Dr. Edmund Halley (1656-1742), the English astronomer, had 
le arned, through a friend, of DeLagny's results for -^08"+^ and 
V a* + b, and this inspired him to investigate approximation meth- 
ods for such expressions. It is amusing to know that a garbled 
version of these results had reached Halley. For the cubic, for in- 
stance, DeLagny found that — |- A/ 'L _|- < \/ ^^ j^ f, <; 
- ^ 4 3a 

a /a? h \ , o-b ^ J/ , ab + a 

--\- 4/-+ ~;a.\soa + —<^a^^b <a+ 

It reached Halley in the form a + < s/ns _i_ 7, < - -|- 

3a^ + b '2 

?!-(-_. As no demonstration accompanied the report of 
4 3a 

DeLagny's work, he derived these statements by a method of his 
own ; this was analogous to that of the French scholar and became 
the basis of his later investigations." He published his results in the 
Philosophical Transactions for' 1694 in an article entitled "A nevv 
exact, and easy method of finding the roots of any equations gener- 
ally, and that without any previous reduction." In this article he 
differs from Raphson and Halley in that he gives frank and generous 
recognition to the work of his contemporaries. 

His approximation for v a^ _)- & is as follows : 

(i) Suppose a' -|- & to be a non-cubic number; also suppose that 
a> + b = {a + ey, approximately, so that a^ < (a + e)^ < {a + i)^. 

Then (a + e)' = a' + ^a^e + sae^ + e' = a^ + b and (a -|- i)' = 
a' + ^a" -|- 30 -|- I ; hence e^< i . 

" Philosophical Transaction of the Royal Society of London, No. 210 (1694), p. 136; 
abridged English translation by Hutton, Shaw, and Pearson, III (1683-1694), pp. 640- 
49. Hereafter referred to as Phil. Trans. (Abr.). 



i. 
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Also b = 30'e + 3ae' + e»; or b = 2,a^e + 3ae^ approximately. 
Since a'e is so much greater than ae', he thinks of e as equal to 

3a^' 

6 6 



From b = 30'e + 3ae^ he gets e = 



a& 



3a' + 3ae . , ('^'\ 



3'^' + 3'^V3a^ 



3a' + 6 



a& 



Therefore Va'+6 = o + e = a+ , . , • This is the rational 

2,a^ + b 

form, which by his handling of e is slightly less than a + e. 

(2) Since b = 3a'e + 3aeS it follows that e' + ae = — ; 

completing the square and finding roots gives the radical form 



a -\- e = — V \/ — -\ . which is a little too large. Similar proc- 

2 '^4 30 

ab , a / n? h 

esses give a < Va^ -b < ' + i/ . The other 

3a* - b 2 y ^ 3a 

results reported by DeLagny are then verified. 

So far Halley has only re-discovered DeLagny's results. But now 
he takes the important step which constitutes his own real contribu- 
tion to facilitating the operation of Newton's method. First he 
derives a general form for all pure powers ; thereupon he extends the 
method of computation as well as its generalization to affected equa- 
tions. It was this important feature that led to Taylor's important 
discovery. Halley says : "While I was considering this matter, I fell 
upon a general method of forms for any power whatever, which 
is elegant enough, and which I cannot prevail on myself to conceal, 
since in the higher powers they at least triple the known figures of 
the root. And these forms, both rational and irrational, proceed as 
follows : 
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</^n> = ^a + ^\aa + A, or a + "^ 



2a^ + -b 
2 



= 5»V[ 



6 , afc 



v^a' + b = — a+ A/ —aa -\ , or a -\- 

4 3<x 3^3 + i, 



3 K 9 



& , ab 



\/a* + b = - a + A/ - aa +—., or a + 



6a 



4a* + h 



=^"+1^ 



b , ab 



</a'+ b = ^a+ A/ — aa + -„ or a + 



6 loa' . , 4, 

50' + ^b 



b , ab 



y/a^ + b = -a + A/ — aa -\ ., or a + 

2 

And so on for the higher powers." 

For the purpose of extending this method to affected equations of 
any degree he builds up what he calls a "general analytical spe- 
culxrni." This is a table of powers to aid in the transforming of the 

equation cz + dz' + fz^ + =0 by letting z = a + e. 

The resulting coefficients of e, e', e', etc., are represented by s, t, u, 
etc. Following is his 

Table of Powers 

s t u w X y 

cz = ca + ce 
dz^ = da^ + 2dae + de^ 

fz' = fa» + 2>fah + 3/ce^ + ffe^ 

gz* = ga* + 4ga^e + 6ga^e' + 4gae' + ge* 

hz<' = ha^ + sha'^e + loto'e' + lohaV + ^hae* + he^ 

k^ = ka^ + 6ka''e + i^ka'^e^ + 2oJfea'e' + i5*aV + 5feae» + ke^ 

The modern reader recognizes the sums of the coefficients in the 
different columns to be the same as in the expansion /(z) = /(a + e) 



= y(a) + y'(a)e + LS^ e + ■^— ^ e» + .... But Halley denotes the 
2! 3! 
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result in the form analogous to that which he uses for pure powers, 

viz., b = se -\- le"^ ■}- ue' + wt* + xe^ + ye* +, where b = f{z) —f(a). 

Operating with this equation he finds that the approximate roots 

sb 

of z are: 2 = a ± , in the rational form, and z = a ^ 

s^ ±tb 



\--A 
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in the radical form. 



Among his numerical illustrations is the equation z^ — 802' + 
19982^ — 149372 + 5000 = o,'* taken from Wallis. Raphson uses 
the same illustration. Even a cursory examination shows the merits 
of Halley's General Analytical Speculum. The arrangement is still 
used in some modern texts on the theory of equations. 



6. Taylor's application of the calculus to 
halley's coefficients 

The next improvement came from the field of the differential 
calculus. Brooks Taylor (1685-1731) had shown in his Methodus 
Incrementorum that 

f{a + e) =f{a)+ f (a) e + ^-^ e' + ^-^ e^+ . . . (Taylor's 

2! 3! 

Series) . 

In an article in the Philosophical Transactions for 171 7" entitled 
"An attempt towards the improvement of the method of approxi- 
mating, in the extraction of the roots of equations in numbers," he 
shows the identity of the coefficients in his series and the coefficients 
in Halley's General Analytical Speculum. Since his series is general, 
applicable to polynomials as well as binomials, to fractional as well 
as integral exponents, to transcendental as well as algebraic expres- 
sions, he observes that Halley's work can easily be extended into 
these fields by aid of the said series. He illustrates this by solving 
the equations (3'^+i)' + 3'— 16 = and log y = 0.29. 

" Phil. Trans. (Abr.), p. 648; Wallis, Algebra, pp. 233-35; Raphson, pp. 25-26. 

" Phil. Trans. (Abr.), VI (1713-1723) pp. 299-304; Phil. Trans., Vol. 30, pp. 610- 

22. 
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7. FURTHER SIMPLIFICATIONS OF NEWTON's METHOD 

In modern text-books we find the Newtonian method of approxi- 
mation stated thus: if /(x) = o, and a is a first approximation, then 

X = a - is a second approximation. Newton and Raphson 

had used these values, but obtained them by transforming the equa- 
tion. Halley used the quadratic part of the transformed equation; 
Newton did the same at times. Brooks Taylor also used the quad- 
ratic, but obtained his coefficient by differentiation. The first 

writer to use the linear differential formula x = a — was Thom- 

as Simpson in his essay, "A new Method for the Solution of Equa- 
tions in Numbers" (1740).^" He states the formula rhetorically and 
illustrates it by this example: f{x) = 3003c: — x' — 1000 = o; x = 
a -f- .4; by trial a is found to lie between 3 and 4; try a — 3.5 ;/'(»:) 
= 300 - 3x^; /(a) ^/'(a) = A; /(3.5) = 7.125; /'(3.5) = 263.25; 
A = 0.027; 3.5 - 0.027 = 3473- 

Proceeding as before, letting a = 3.473, we have x = 3.47296351. 

Simpson uses Newton's method in solving simultaneous equations. 
His fifth example is a set of exponential equations: x'' + y — 1000 
= o and x^ + y^ — 100 = o. 

De Courtivron (1744), Euler (1755), and Lambert (c. 1770) also 
applied the calculus to the solution of higher equations.'" In his 
algebra ^^ (1770), Euler discards the higher powers in both dividend 
and divisor during the substitution instead of carrying them along 
until the division takes place. Thus in solving the equation x'+ 
ax^ -\- bx -\- c = o, he sets x = n — p, x^ = n^ — 2np, x' = w' — 

, 2>, f tu { ■ 1, c ^ A n^ + an^ + bn + c 
^n^p; from these four expressions he nnds p = 

3«^ -1- 2an + b 

Edward Waring (1770) extended the Newtonian approximation 
process into the realm of complex roots. The functional expression 

X = a — ^-^ is due to Lagrange ^^ (1798). 
fia) 

2" Simpson, Essays on Mathematics, London, 1740, p. 81. 

21 Euler, Institutiones calculi differentialis, sections 224, 234, 235; Histoire de 
I'Acadlmie, 1744, p. 405-14. 

^ Euler, Algebra, I (1770), Chap. 16. 

■' Lagrange, Traili de la Rlsolulion des Equations Niimeriqtus (1798), Note XI. 
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8. SUMMARY 

1. The method of solving equations of higher degree by suc- 
cessive substitutions in dc ived equations was originated by Newton 
(1669) incidental to his getting integral expressions for his work in 
areas. He performed a transformation for every new supplement to 
the root. His method was first made public in the algebra of Wallis 
(1685). 

2. Raphson (1690) worked Newton's method into a system. 
Following the example of Harriot he derived canonical forms to 
simplify calculation. In his plan only one transformation and one 
divisor are needed for each equation. 

3. DeLagny (1691) derived notable approximations in the solu- 
tion of pure equations by a special adaptation of Newton's method, 
and worked out an elaborate set of tables for the solution of numer- 
ical equations. 

4. Halley (1694) generalized DeLagny's method for pure equa- 
tions and extended it to affected equations ; in his table of powers he 
gave a general formula for the coefficients of transformed equations. 
This made obsolete Raphson 's canonical forms and DeLagny's 
tables. 

5. Taylor (1717) discovered that Halley's coefficients were the 
same as the coefficients in Taylor's series. This made possible the 
direct use of the calculus in finding roots. 

6. Our present standard linear approximation formula x = a — 

was first used by Simpson (1740); it received its conventional 

form through the writings of Lagrange (1798). 



VIII 

CERTAIN EFFECTIVE BUT NON- 
PRACTICABLE METHODS 

I. THE ASCENDANCY OF NEWTON 's METHOD 

Newton's method of approximating roots by substituting values 
in equations of an infinite number of terms did not immediately 
upon its publication in 1685 supersede Vieta's method of finding 
roots by evolving the successive digits. The latter was still used by 
Dech&les in 1690 and appears as late as 1702 in a new edition of 
Oughtred's Clavis. It seems as if the transition from the one to the 
other must have been facilitated by the eclectic character of Raph- 
son's arrangement. Halley's general formula increased the popular- 
ity of the new method. Its ascendancy is clearly shown in an article 
written by Wallis in 1694. After Taylor's discovery (1717) it had no 
serious rival old or new, until the publication of Horner's method in 
1819. But new methods sprang up to challenge it in this period, and 
we shall describe the three that gave the most promise. 

2. ROLLE's METHOD OF CASCADES 

A considerable contribution to the location of roots was made by 
Michel Rolle in his two works Traite d'Alghhre (1690) and DSmons- 
tratio d'une MSthode pour Resoudre les Egalitez de toutes Degrez 
(1691). He demonstrated that there cannot be more than one real 
root of f{x) = o between two successive real roots of /'(x) = o 
(Rolle's Theorem). He invented a method of finding roots, called 
the "Method of Cascades"; it consisted in locating the roots by first 
finding the limits ("hypotheses") of the roots of the different cas- 
cades. Between the roots of each cascade lie the roots of the next 
higher cascade, and so on.' It was an anticipation of the method of 
using successive derivatives (Rolle's cascades) now used in calculus 
for locating roots. 

With Rolle an equation is "piepared" when it has the form 
x" — Uix"'^ + aix"'^ — flax""^ + =0 where the a's are posi- 

' Traiti d'Algebre, p. 124, p. 103; p. 118; Dlmonslratio, article VI (New York Public 
Library); Cajori in Bibl. Math. XI (1910-11), p. 300-313. 
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tive integers. Then the upper limit of the root is found by dividing 

the largest negative coefficient by the highest exponent, and adding 

unity and as much more as necessary to make the result an integer. 

His Method of Cascades he illustrates thus:^ 

Let x' — 57^^ + 936*: — 3780 = o (third cascade). 

Multiply the terms of the cascade by parts of the piogression 3, 2, 

I, o, respectively; divide by x. Then 3*^ — 114a: + 936 = o (second 

cascade) . 

Multiply by 2, i, o; divide by 2X. Then 3a: — 57 = o (first 

cascade) . 

The "hypotheses" for the second cascade are o, 19 and 39 (or 

'-^ + I). 
3 

The "hypotheses" for the first cascade are o and iq. 

Between o and 19, the second cascade is found to have 12 as a root; 

and another root, 21, between 19 and 39. 

Finally for the third cascade the extreme "hypotheses" are o and 3781 

( or — h I ) ; o, 12, 26, and 3781 are its four "hypotheses." 

The roots of the third cascade (that is, the given equation) can 
then be found to be: 6 between o and 12; 21 between 12 and 26; 30 
between 26 and 3781. 

3. THE METHOD OF RECURRING SERIES 

De Moivre (1667-1754) first conceived of series in which each co- 
efficient bears a given relation to certain of the preceding coeffi- 
cients; such series he called lecurring series (1720).' Daniel Ber- 
noulli (1700-1782) first used them for approximating roots (1728). 
Following is his process : * 

Take such an equation as i = ax + bx^ + cx^ + dx*. Select 
arbitrarily the four terms A, B, C, D. Then E, F, G, H, etc., con- 
stitute a recurring series if 

' Traiti d'Algibre, pp. 127-28. 

'Phil. Trans., Vol. 32 (1722), pp. 162-78; DeMoivre, Miscellanea Analytica, 
London, 1730, Bk. I, II. 

* Cantor. Ill (1898), pp. 621-22. 
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E = aD + bC + cB + dA, 

F = aE + bD + cC + dB, 

G = aF + bE + cD + dC, 

H = aG + bF + cE + aD; 

F F 
and the approximate values of x are—,-, etc. 

F G 

In solving the equation i = — 2x+5x;' — 4:c^ + x*, Bernoulli sets 
A=B = C = D=i, and obtains the series 
0,2, -7,25, -93-341- -1254- 

This gives -,— -,— — , as the increasingly close 

2 7 25 1254 

approximations of the root. The last value, when substituted in 
the given equation gives i = 0.999487, a very small error. 

Bernoulli gave no analytical proof for the validity of his process. 
That was supplied by Euler (1748) and later by Lagrange (1798).' 

A considerable improvement was made by Euler in his algebra 
(1770). He shows the general relation that must exist for the dif- 
ferent powers, and thus immediately gets a recursion formula for 
any equation under consideration.* The recursion formula performs 
a function similar to that of Raphson's canonical form or of Halley's 
general formula. 

Suppose our series p, q, r, s, t, etc., is such that -, -, etc., express 



Of . , T 

approximate values of x. Then: = x;- — x; multiplying,- = x^; 



P 1 
lultip 

g ■ p 



S S t 

also, - = x; hence- = x^; similarly - = oc^, etc. Not only do we 
r p p 

have a substitution table for powers, but we have one with the first 
term as the common denominator. 

As an illustration, let us substitute these values in the equation 
x^ = x"^ -\- 2x -\- I. This gives s = r -{■ 2q -\- p. Using this as a re- 

' Euler's IntroducHo in Analysin Infinilorum, Chap. 17; Lagrange, Traitl de la 
Resolution des Aquations Numlriques, Note VI. 
' Euler's Algebra, Vol. I, Chap. 16. 
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cursion formula for 5 in terms of r, q, and p, we can arbitrarily select 
these to be o, o, i. We then get as our series o, o, i, i, 3, 6, 13, 28 

60, 129, etc. Testing ior x = — in the original equation we get an 

28 

I'x . 129 

error of — - . A closer result would be given by ac = — . 
343 60 

Not all equations admit of the application of the method of 
recurring series. If no common limit results, the equation must first 
be transformed. This happens particularly if the second term is 
lacking. Thus foi the equation a;^ = 2 we get the recurring series 
I, I, 2, 2, 4, 4, 8, 8, etc. A similar situation arises for x' = 2. To 
obviate this, we let a; = y — i. In the resulting equation y' = 
3y^ — 3>' + I we get the recursion formula s — y — sq + p, and 
this gives us the series 0,0, i, 3, 6, 12, 27, 63, 144, 324. This gives us 

y = 5-^ and x = -. Checking will reveal an error of — . 
144 4 64 

4. THE METHOD OF CONTINUED FRACTIONS 

Newton's method had one serious defect: it was inoperative if 
two roots were close together, for then the series would not converge, 
and might after a while actually diverge. Besides this, Lagrange 
pointed out two other defects: the uncertainty as to the exactness 
of each fresh correction, and the failure to give a commensurable 
root in finite terms. Lagrange announced (1767) a new method 
which was free from these defects.' It consisted of three parts: (i) 
a scientific method of finding the integral part of the root instead of 
the empirical methods hitherto employed; (2) a rule for separating 
the roots; and (3) a technique of approximation by the use of con- 
tinued fractions. This last may be explained as follows: 

Suppose that in the equation /(*) = o we have found that there is 
only one root between the integers p and p -\- 1; for convenience let 

it be real and positive. Then x = p + -; and y is greater than 

y 

unity. Substituting the value of x 'mf{x) gives an equation F{y) = 0. 

' Lagrange. Traili, Chap. I, III, IV and Notes IV, V. 
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Since y > i, find its integral value by the prescribed method. 

Suppose it lies between the integers q and g + i. Then y = q + -; 

z 

and so on. Then x = p -\- 



r + .. 



Lagrange illustrates his method by numerical examples. One of 
these is Newton's equation x' — 2x — 5 = o.^ By his methods of 
limits and separation he finds one root between 2 and 3. Let x = 

2 +-; this gives the equation 31' — lo^ — 6y — i =0. Again 

y 

applying the method of limits and separation he finds that y lies 

between 10 and 11. Similarly for y = 10 + - he derives the equa- 

z 

tion 6iz' — 94^^ — 20z —1=0; here i < z < 2. 

Letting z = 1 + - gives 54M' + 25u^ — 8gu — 61 =0: here 

u 

I < M < 2. 

A continuation of this process gives the series 2, 10, i, i, 2, i, 3, 
I, I, 12. 

Hence x = 2 -\- i 



10 + I 



I + I 



I + I 



2 + 



Evaluating this expression by known methods gives the approxi- 
matins series- 22i2344HiiL557673ii307 i64£5 ^^ 

ludung Scries. — , — ; — i — i t r » » , » » cll., 

I ID II 21 53 75 275 349 624 7837 

» Train (1798), Chap. IV, p. 29. 
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That is, X = — ^^—^ = 2.09455149; and the error is less than 
7837 

T 

= 0.0000000163. 



(7837) = 



Though none of these methods gained popular usage, they gave 
rise to by-products that were incorporated by algebraists into 
Newton's method. Among these were Rolle's theorem for locating 
roots and the method of Lagrange for finding the integral parts of 
the root, as well as his rule for separating roots. Waring had 
already (1762) formulated such a rule but it was very little known 
until Lagrange published his method of squared differences. 

5. SUMMARY 

1. Between 1685 and 1819, three potential rivals of Newton's 
method sprang into being: Rolle's Method of Cascades (1690), 
Bernoulli's (1728) and Euler's (1748, 1770) method of recurring 
series, and Lagrange's method of continued fractions (1767). 

2. None of these methods, though scientific and effective, became 
generally adopted, and Newton's method was supreme from about 
1694 to 1819. 



IX 

HORNER'S METHOD; SIMILAR METHODS BY RUFFINI 

AND BY THE CHINESE OF THE 

THIRTEENTH CENTURY 

I . A GENERAL VIEW 

No method, not even that of Vieta, had been so universally 
adopted as that of Newton. Of its two latest rivals the method of 
recurring series, though as presented in Euler's Algebra not par- 
ticularly laborious, never had its possibilities exhausted; and La- 
grange's method of continued fractions, though effective and perfect 
in theory, was too laborious for popular adoption. At the beginning 
of the 19th century, therefore, Newton's method held undisputed 
sway. Its one drawback had been its insecurity. But the success of 
Waring and Lagrange in effecting methods for the separation of the 
roots had to a large degree eliminated this defect. The discoveries 
of Fourier and Budan, early in the nineteenth century, whereby the 
number and nature of the roots could be detected for any interval, 
gave Newton's method a much more scientific basis. ^ The mechan- 
ism of its operation had been so greatly improved by Halley, Taylor, 
and Simpson, that in this respect it has to this day remained without 
a peer. 

It has had only one serious rival. It is that form of approxima- 
tion by the evolution of digits made known to the Western world by 
William Horner in 1819. Like the methods of Vieta and Newton, it is 
also a process of exhaustion. Horner's process resembles that of 
ordinary division or of current methods of finding square and cube 
roots. That element was not new, having been the chief principle in 
Vieta's method. But Horner's way of doing this by the gradual 
diminution of the roots in a succession of transformed equations, 
was new to Europe. To the world it was old, however; it had been 
known and used in at least three mathematical centers before 1819, 

' As early as 1768 J. Raytn. Mouraille had successfully formulated a plan for insur- 
ing the security of Newton's method; but his work was unknown at the time, and has 
only recently been noticed. See article by F. Cajori in Bibl. Math. (3) XI (191 1). pp. 
132-37- 
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— China, Japan, and Italy. In this chapter we shall give a short 
survey of the different forms this method has assumed. 



2. THE METHOD USED IN CHINA AND JAPAN 

As early as the first half of the 7th century Wang Hs'iao worked 
with numerical cubic equations in his Arithmetic in Nine Sections? 
He found the roots by using a method similar to that of finding the 
cube root of numbers. He attempted no equations of degree higher 
than the third. 

Not till six centuries later do we find any work with equations 
higher than cubics.^ Then we find three men, living almost at the 
same time, systematizing a method for the orderly evolving of the 
digits of the root. The first, Ch'in Chiu-shao, in his Nine Sections 
of Mathematics (1247), solved equations of the 6th, 7th, 8th, and 
even higher degrees. He worked on tables, with wooden rods to 
represent his coefficients, and arranged his work in columns as did 
the Hindus for roots of numbers. In writing his equations with the 
right hand member equal to zero he anticipated Harriot by over 
three hundred years. Li Yeh, the second pioneer in this work, in 
The Mirror of the Mensuration of Circles (1259), explained the 
process more fully; Ch'in had chiefly given rules. The third member 
of the trio was Chu-chieh. His form of the process, as presented in 
his Introduction to Mathematical Studies (1299), found its way into 
Japan in the seventeenth century. The names given to this method 
are quite suggestive. Ch'in Chiu-shao (1247) called it "Harmoni- 
ously alternating evolution"; Yang Hwuy, another writer of the 
same period, called it "Accumulating involution." 

To show the steps in this process, invented in China in the 
thirteenth century and introduced into Japan in the seventeenth 
century, we shall give the solution of a problem in which the root is 
12, taken from the eighth book of Sato Moshun's Tengen Shinan 
(1698). The arrangement ^ is Occidental. The solution is as follows: 

' Yoshio Mikami, The Development of Mathematics in China and Japan, New Yo rk 
1914; pp. 53-56. 

' D. E. Smith and Yoshio Mikami, A History of Japanese Mathematics, Chicago, 
1914; pp. 49-56; Yoshio Mikami, pp. 63-89. 

* Given by D. E. Smith and Y. Mikami. The Chinese arrangement is given in full in 
Smith and Mikami's Hist, of Jap. Math., pp. S3-55. 
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Given x* + 4a:' — 236^;=' — 432a; + 11520 = o 
1+4 - 236 - 432 + 1 1520 
10 140 — 960 — 13920 



I 


14 

10 


- 96 

240 


- 1392 
1440 


— 2400 


I 


24 

10 


144 
340 


48 


— 2400 
2400 


I 


34 
10 


484 


48 
1 152 





I 


44 
2 


484 
92 


1200 
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Among Japanese scholars to write on the Chinese method ' 
of approximation were Sato Seiko (1666), Seki Kowa (1642-1708), 
Sato Moshun (1698), Sakabe (1803), and Kawai (1803).' 

3. THE METHOD INVENTED BY RUFFINI 

In 1804 a solution very similar to the approximation process of 
the Chinese was effected by Paolo Rufifini.^ By means of the differ- 
ential calculus he worked out a theory for transforming one equation 
into another with roots diminished by a certain constant. In his 
revisions of 1807 and 1813 he used algebraic processes for the same 
purpose. His arrangement differs from that of Horner in that the 
coefficients of the transformed equation appear in the extreme right 
hand column, whereas in Horner's arrangement they appear in a 
diagonal line. A single solution will tend to show the similarity of 
Ruffini's work to Horner's process and to that of the Chinese: 

' Whether the process used by the Chinese was actually the same as that used by 
Horner is questioned by Professor Gino Loria, of the University of Genoa, Italy, in an 
article published in the BoUetino della Matkesis for April, 1920, and abridged by Pro- 
fessor R. B. McClenon, of Grinnell College, Grinnell, Iowa, in The Scientific Monthly 
for June, 1921. In discussing Ch'in Chin-Shao's solution of the equation — «* + 763, 
200*'— 40, 642, 560, 000 = o he says: "We still have too meager information as to 
the details of the work for us to be able to affirm confidently that Horner's method was 
known to the Chinese in the 13th century; we can only say that this method, or one 
practically identical with it, was known at the time, and we must await further evidence 
before affirming or denying the priority of the Chinese in its discovery." 

' Smith and Mikami, pp. 53, 86, 116, 212-15. 

' F.. Cajori, "Horner's Method of Approximation Anticipated by Ruffini," in Bull. 
American Math. Soc, XVII (1911), pp. 409-414. 
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"Given the equation ^x^ — 6x^ + 3a:' — ^x' — ^ + 8 = o; re- 
quired, to transform it into another in which y = x — 6. Operate 
in the manner here set forth : 

4- - 6, 3. - 5. - 4. 8, 

4, 18, III, 661, 3962, 23780 

4, 42, 363, 2839, 20996 

4- 66, 759, 7393 

4, 90, 1299 

4, 114 

4 
and we obtain the transformation 

4^5 + 114/ + i299j'S + 73933-2 _|_ 209963) + 23780 = o." 

For Newton's classic example x' — 2a; — 5 = o he finds the root 
to 8 decimal places, by this method. 

Ruffini's method remained unnoticed by his own people. It was 
the story of Vieta repeated. But there was no Harriot or Oughtred 
to bring the message to England this time. For Britain had her own 
discoverer in this field. It is a coincidence that every successful 
approximation method for numerical equations has been either in- 
vented or sponsored in Great Britain. 

4. Horner's method 

The theory and mechanism of Horner's method are so well known 
that we need not discuss them here.' We shall review only the out- 
standing points in its history. 

It is safe to say that neither Horner nor Ruffini knew of the other's 
work, nor was either acquainted with the methods of the Orient. 
Horner's paper was read before the Royal Society of London on July 
I, 1819, and was published in the Philosophical Transactions for that 
year.' 

In developing his theory, he used Taylor's Theorem and Arbo- 
gast's derivatives (Ruffini had used ordinary derivatives), and 

* A thorough discussion of its analytic and operational phases is to be found in J. R. 
Young's Theory and Solution of Algebraical Equations of Higher Order, London, 1843. 
Chapter XII. A vigorous challenge to rival English claimants for priority over Hor- 
ner is to be found in De Morgan's article already referred to. 

' W. H. Horner, "A new method of solving numerical equations of all orders by 
continuous approximation," in Phil. Trans., Vol. 109 (1819), pp. 308-335. 
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arranged his results in a "General Synopsis" in much the same way 
as Halley built up his General Analytical Speculum. He obtained 
his transformations by using the calculus, but says that algebraic 
processes can be used for this purpose. Like Ruffini he used ordi- 
nary algebra in his revisions of 1830 and 1845. 

His first illustration is the classic example of Newton. The work 
is arranged differently from that which we see in text-books on al- 
gebra. His schema resembles Wallis's arrangement of Vieta's process 
and our present arrangement in finding the cube root; he even 
divides the resolvend into periods. 

Horner's method met with quick response in England, due in part 
to its enthusiastic advocates, J. R. Young and A. De Morgan. It is 
today a strong rival to Newton's method in both England and the 
United States. It is widely used, though less than Newton's method, 
in Germany, Austria, and Italy. In Fiance, Newton's method has 
held undisputed sway. Horner's method, when combined with the 
theorems of Budan and Fourier and of Sturm, is effective, secure, and 
expeditious. 

5. SUMMARY 

1. In the 13th century Ch'in Chiu-shao (1247), Li Yeh (1259), 
and Chu-chieh (1299) employed a method of approximation vir- 
tually the same as Horner's method. In the 17th century it was in- 
troduced into Japan. 

2. Paolo Ruffini (1804) invented a similar method in Italy, which 
was soon forgotten. 

3. Horner published his method in 1819, and it soon became 
widely used in England and later in the United States, and to a less 
degree in Germany, Austria, and Italy. It was never adopted in 
France. 



X 

GENERAL SUMMARY 

1. The earliest known method of finding the roots of numbers 
and of solving formal equations is that of Single False Position. 
This method was used by the Egyptians, the Greeks, and the Hindus. 

2. The method of averages was probably used by Archimedes; 
a definite exposition of it is given by Heron of Alexandria. Its appli- 
cability was restricted; but it was often incorporated into other 
methods, as in Chuquet's Rule of Mean Numbers. 

3. Methods of Double False Position for finding the roots of 
numbers seem to have been used by Archimides; a specific process 
was described by Heron of Alexandria. In algorithmic work double 
false position was used by the Arabs and Leonardo of Pisa. By the 
Renaissance writers these processes were extended to the solutions 
of numerical higher equations. Though eclipsed by the methods of 
exhaustion invented by Vieta, Newton, and Horner, methods of 
double false position have never become obsolete. Two main plans 
for finding mean values have been followed: (i) the averaging or 
analogous devices, as illustrated by Chuquet and DeLagny; (2) in- 
terpolation by proportional parts as used by Leonardo of Pisa, al- 
Karchl, Burgi, and Pitiscus. 

4. Methods of exhaustion have been used the most extensively 
of all approximation methods. The earliest method of exhaustion 
recorded is that used by Theon of Alexandria for irrational square 
roots by reversing the binomial formula. By a special adaptation of 
this method Aryabhata, Brahmagupta, Sridhara, and Bhaskara 
worked out a process for evolving the roots of numbers digit by digit 
(Evolution). This was adopted by the Arabs and by them trans- 
mitted to Christian Europe through the writings of Leonardo of Pisa 
and of Sacrobosco. Oui present standardized process of finding roots 
of numbers is largely the work of Peurbach and Cardan. 

Vieta extended this method to numerical higher equations. His 
method was popularized and made more practicable by Harriot, 
Oughtred, and Wallis, and remained the prevailing method until 
superseded by Newton's method. 

5. Newton's method of exhaustion, based on equations whose 
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terms constitute an infinite series, was conceived as an aid to his 
work in the calculus. The place of this method as a constituent part 
of algebra is largely due to its modification and systematic develop- 
ment by Raphson. Its technique was improved by DeLagny, 
Halley, Taylor, and Simpson and its scientific basis was clarified and 
strengthened by Lagrange, Fourier, Budan, and Sturm. No other 
method of approximation has come up to it in general popularity. 

6. A third method of exhaustion, employing continued fractions, 
was invented by Lagrange. It is a theoretically perfect method, but 
too laborious for practical use. It has retained chiefly an academic 
interest. 

7. Rolle's Method of Cascades and the method of recurring 
series developed by Bernoulli and Euler have been used only slightly. 
The methods invented by Collins, Fontaine, and some of those used 
by DeLagny were never taken up into the mathematical activities 
of the world. 

8. A fourth method of exhaustion in solving numerical higherequa- 
tions is both the oldest and the youngest, — oldest, because it was 
used by the Chinese in the thirteenth century, youngest, because it 
was re-invented and given to the Western world by Ruffini and 
Horner in the twentieth century. It is a process of evolution of 
digits resembling that of Vieta; however, the digits are evolved, not 
by the inversion of a binomial formula, but by the transformation of 
equations. Horner's method has made progress in Italy, Germany, 
and Austria, but is hardly ever used in France; in England and the 
United States it is a close lival of Newton's method. 



APPENDIX 

Note i 

The solution of numerical higher equations by means of graphs will often 
give as near an approximation as is desired. ' Where greater accuracy is 
desired a combination of the graphic method with one of the standard methods 
of exhaustion or of double false position is generally advantageous. The 
application of the calculus, such as is ordinarily used in curve tracing, facili- 
tates the work of approximation; in isolating the roots, for example, one can 
often use Rolle's thiK)rem instead of the more elaborate theorems of Fourier, 
Budan, and Sturm. 

Whether tables will ever be used in the solution of numerical equations is 
for the future to tell. DeLagny advocated their use and worked out a 
comprehensive set of tables similar to our modern integration tables in his 
Analyse Generate; but, as his voluminous work has been very little known, his 
idea has never been tried out. 

Note 2 

Teachers of algebra often fail to recognize how extensively the methods of 
double false position are used in mathematical thinking. Notice how similar 
to Btirgi's method of solving equations is the following common method of 
determining the autumnal equinox: As given in the Almanac of 1899, the 
sun's declination at approximate noon (Greenwich), Sept. 22, was + 0° 18' 
8" . 7 and at noon, Sept. 23, it was — 0° 5' 13" . 9. The difference in declina- 
tion being 23' 22" . 6, the equinox is found, by interpolation, to be 18' 18" . 7 
-i- 23' 22" . 6 or 0.7762 of a day after noon, Sept. 22. Take also the Law of 
Charles in physics, by which we know that if the temperature of a gas be 

lowered 1° C. its volume will be decreased by — of its volume at o* C, and if 

273 

its temperature is lowered 10° C. its volume will decrease by — of its volume 

273 

at 0° C. From these two positions the volume at any temperature, say 
— 100° C, for which the gas does not become a liquid, is calculated. The 
Fahrenheit-Centigrade thermometer problems, the interpolations for trigo- 
nometric functions, and many problems in proportion are more familiar illus- 
trations. 

' Consult Arthur Scliultre. Graphic Algebra. New York, 1908. 
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The formal and informal interpolation processes used in astronomy, physics, 
trigonometry, and other computational fields outside and inside the class- 
room constitute a big share of the actual calculation. There seems to be no 
valid reason, then, why the student of algebra should not have his attention 
called to the use of such methods in the solution of numerical higher equa- 
tions. While less refined than methods of exhaustion, they are very usable for 
securing results correct to one or two decimal places. An interesting text- 
book exposition is that given by Professor Isaac Dalby in the second volume 
of his A Course of Mathematics, designed for the use of the Officers and Cadets of 
the Royal Military College (London, 1806), pp. 175-185. For good illustrations 
of the use of these methods in recent texts, see Wells's University Algebra 
(Boston, 1884), pp. 417-19, and F. L. Griffin's Introduction to Mathematical 
Analysis (New York, 1921), pp. 37-38. 

To find the mean between the two false f)ositions we generally use one of 
the following three ways: (i) finding the arithmetic average, or using some 
such special formula as Chuquet's Rule of Mean Numbers; (2) interpolating 
directly by proportional parts; (3) using the interpolation formula 

* = .where oi and 02 are the approximate values and ei and ej 

ei — 62 

are the corresponding errors. 

Note 3 

The solutions of numerical higher equations first became a mathematical 
need in the work of geometry, trigonometry, and astronomy. However, they 
enter into expert practical calculation more frequently than the average 
student realizes, and the author feels convinced that the insertion of such 
problems would vitalize the methods of approximation given in our college 
algebras. The following representative problems from the fields of physics, 
mechanics, economics, and engineering, selected from the more progressive 
text-books, will serve to illustrate:''' 

(i> The real root of the equation x* — 383c — loi = o gives the diameter 
of a water pipe 200 feet long, discharging 100 cubic feet per second under a 
head of 10 feet. 

Result: diameter is 2.92 ft. 

(2) The equation ifi — 19200W -|- 21 1200 = o gives the speed in feet per 
second on a i-inch manila rope transmitting 4 horse-power, under a tension 
of 300 pounds on the tight side. 

2 Merriman and Woodward, Higher Mathematics, p. 13; Rietz and Crathorne, Col- 
lege Algebra, p. 140; Karpinski, Benedict, and Calhoun, Unified Mathematics, p. 404; 
F. L. Griffin. Introduction to Mathematical Analysis, pp. 37-42. 337-342, More prob- 
lems of a similar nature may be found in these same references. 
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Result: speed is 11.07 feet. 

(3) From the American Report on Wholesale Prices, Wages, and Trans- 
portation, for 1891, the median wage in dollars is given by X of a value of x 

in the equation 2561 - = 6972 — + 657 — x — -ii, — x^+ — x' — x*. 

2 128 18 64 48 128 

Result: wages were $1,536. 

(4) In Archimedes's problem of the section of a sphere, if k represents the 
ratio of the larger to the smaller segment, the distance of the plane section 
from the center of a sphere 10 inches in diameter is given by the equation 



30OX + 2000 



(rrl) = " 



(5) A sphere of yellow pine i foot in diameter floating in water sinks to a 
depth X given by the equation 2x' — ^x^ = 0.657. 

Result: depth is 0.606 foot. 

(6) The equation d* + 320^^ — 340^ — 4290 = o gives the smallest safe 
diameter (in inches) for a certain steel shaft. 

Result: the diameter is 4. 11 inches. 

(7) The equation 2x^ — gx' + 24 = o gives the depth of immersion in 

Q 

water of a sphere 3 feet in diameter and with a specific gravity -. 

Result: depth is 2.38 feet. 

(8) The equation x* — iiSx'' + i6ox — 115 = o gives the length of the 
longest rectangular panel I foot wide which can be fitted diagonally across a 
door 4 feet wide and 10 feet long. 

Result: length is 10.17 feet. 

(9) The "index of correlation" between the eye-colors of a certain group of 
people and of their great-grandparents is approximately the root of the equa- 
tion 0.024X* -|- 0.137X' -f- 0.035^^ + X — 0.225 = o which lies between 
and I. [C. B. Davenport] 

Result: index is 0.22. 

(10) Using Jupiter's radius, 45090 miles, as the unit of distance, the equa- 
tion x' — 5x' + 6.27396X — 0.060385 = o gives the greatest and least dis- 
tances of Jupiter's Fifth Satellite from the center of the planet. 

Results: distances are 2.48 and 2.51 units. 

(11) A magnet placed with its ends in a ''magnetic meridian" will neutralize 
the earth's magnetism at certain points. To calculate the position of these 

, , . 2000OX 

pomts m a certam case, it was necessary to solve the equation — = 0.2. 

{x^ — 100)' 

Result: x = 47.82 and 0.09. 
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Note 4 

Bearing upon source material relating to the solution of numerical equa- 
tions, the reader may be interested to know that in New York libraries are 
found certain works very rare in this country. Among these may be men- 
tioned the manuscript copies of the Lilawati, the Kholasat-al-Hisdb, and of 
Chuquet's Triparty in the private collection of Professor David Eugene 
Smith. In the Columbia University library there is a manuscript on algebra 
by Scheubel, in his own handwriting, prepared for the press but never pub- 
lished; also a copy of Vieta's collected works published by van Schooten in 
Leyden in 1646, and a copy of J. Hume's exposition of Vieta's algebra, 
published in Paris in 1636. The latter work seems to have been unknown to 
the mathematical world, as we find no references to it in any published history 
of mathematics.' Hume anticipated the "Newtonian divisor" by an analogous 
device, but his work was evidently unknown to his contemporaries. In the 
New York Public Library are found three comparatively rare books: Hdri- 
gones' Cursus Mathematicus (Paris, 1642), which popularized Vieta's method 
in France; Michel Rolle's Demonstratio d'une MSthode; and Raphson's 
Analysis Aequationum Universalis, edition of 1697, bound with his almost 
wholly unknown metaphysical treatise, De Spatio Reali. 

' It is referred to in Professor D. S. Smith's forthcoming history of mathematics. 



VITA 

Martin Andrew Nordgaard was born January 29, 1882, near 
Northwood, Iowa. He graduated from Northwood High School in 
1898, and in 1903 he received his Bachelor of Arts degree from St. 
Olaf College, Northfield, Minnesota. From 1903 to 1909 he served 
as teacher and administrator in the secondary schools of Wisconsin 
and Minnesota, with the exception of one year spent in the study 
of education at the University of Minnesota. He was the head of the 
department of mathematics at Columbia Lutheran College in 
Everett, Washington, from 1909 to 1912. The following year he 
was a graduate student of mathematics at the University of Chicago. 
The three years from 1913 to 1916 he was instructor of mathematics 
at the University of Maine. From this institution he received the 
degree of Master of Arts in 1914. During the year 1916-17 he served 
as interim head of the department of education at St. Olaf College. 
In 1917-18 he held the Harrison scholarship in mathematics in the 
graduate school of the University of Pennsylvania. In the fall of 
191 8 he became a member of the mathematics staff of Grinnell Col- 
lege, Grinnell, Iowa. In February, 1921, he was granted leave of 
absence for the purpose of further study. Since then he has been a 
graduate student in Teachers College, Columbia University, major- 
ing in the department of mathematics. During the current year he 
has held the position of lecturer in mathematics at Columbia Uni- 
versity. 



'^m: 



'•'Ti*^ 












r^^ 



